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1.0  INTRODUCTION 


The  purpose  of  this  research  is  to  study  a  nondimeusional  measure  of  struc¬ 
ture-control  system  perrormance  which  has  the  potential  to  characterize  both 
quantitatively  and  qualitatively  the  designer's  ability  in  dealing  with  so.ie  ot 
the  problem  areas  such  as  assessment  of  control  spillover  effects,  model/con¬ 
troller  order  reduction,  input  configuration,  and  the  interaction  between  struc¬ 
tural  and  control  variables  from  the  structure-control  system  point  of  view.  The 
nondimonsional  measure  is  defined  as  the  efficiency  of  the  system. 

The  concept  of  efficiencv  is  widely  encourtered  in  thermal  and  thermo¬ 
mechanical  sciences.  In  these  fields  an  efficiency  is  defined  as  a  nondimen- 
sional  ratio  of  two  scalars  which  represent  energy  or  energy  related  quantities. 
Typically,  one  of  the  scalars  characterizes  a  theoretically  ideal,  but  physically 
irrealizable  process,  and  the  other  characterizes  an  actual  physically  realized 
process.  The  difference  between  these  two  scalars  represents  waste  of  the  total 
available  energy  or  energy  related  quantity  in  realizing  the  actual  physical 
process  and  is  regarded  as  an  irreversibility  inherent  in  the  physical  system. 
Many  different  forms  of  efficiency  are  defined  ir  thermomechanical  sciences 
depending  on  how  the  theoretically  ideal  and  the  physically  realized  actual 
processes  are  identified  in  a  certain  application.  Some  examples  are  propulsion 
efficiency,  heat  engine  and  heat  pump  efficiency,  adiabatic  compressor 
efficiency,  brake  efficiency,  overall  efficiency,  etc.  (Ref.  9).  In  thermo¬ 
mechanical  disciplines  an  important  objective  is  to  design  a  system  with  high  or 
maximum  efficiency  consistent  with  the  physical  constraints.  In  practical 
engineering  terms,  the  thermomechanical  systems  which  use  (or  convert)  the 
largest  fraction  of  the  total  available  energy  in  realizing  the  physical  function 
of  that  system,  thus  causing  the  minimum  waste  of  resources,  are  desired. 
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The  utility  of  efficiency  concepts  as  analysis  and  design  tools  in  thermo- 
mechanical  systems  is  well  established  and  provides  valuable  physical  insight  to 
the  working  of  the  system.  The  question  arises  whether  such  a  time  tested 
concept  in  thermomechanical  discipline  can  be  extended  to  the  field  of  distrib¬ 
uted  parameter  structure-control  systems  and  yield  comparable  practical  value  and 
physical  insight  for  the  analysis  and  design  of  such  systems.  This  investigation 
is  a  quest  in  that  direction  and  presents  a  conceptual  framework  to  establish  the 
usefulness  of  efficiency  concept  for  structure-control  systems. 

For  structure-control  systems,  the  efficiency  concept  is  defined  as  the 
ratio  of  two  control  functionals  pertaining  to  the  particular  structure-control 
design.  The  control  functionals  are  judiciously  defined  to  represent  the  average 
control  powers  consumed  during  the  control  period.  For  the  purpose  of  defining 
efficiencies,  four  quadratic  control  power  functionals  are  relevant.  Out  of 
these  four  control  powers  a  relative  model  efficiency  and  a  global  efficiency  can 
be  defined  for  a  structure -control  system  with  maximum  possible  percent  effi¬ 
ciency  of  100  for  each.  The  two  efficiencies  are  related  by  a  modal  efficiency 
coefficient. 

Relative  model  efficiency  is  defined  as  the  ratio  of  a  quadratic  modal  con¬ 
trol  power  functional  to  a  quadratic  real  control  power  functional  associated 
with  any  control  design  model.  Physically,  relative  model  efficiency  represents 
the  fraction  of  real  control  power  expended  on  the  distributed-parameter  system 
(DPS)  usefully  absorbed  by  the  control  design  model.  The  relative  model 
efficiency  can  be  defined  for  control  of  the  DPS  by  both  spatially  discontinuous 
(discrete)  and  spatially  continuous  input  configurations.  If  the  dimension  of 
the  control  design  model  and  tl  DPS  system  are  identical  then  the  relative  model 
efficiency  by  definition  becomes  100%.  Consequently,  the  relative  model 
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efficiency  is  predominantly  an  indicator  of  the  effect  of  finite 
dimensionalization  of  the  control  design  model  in  using  the  available  control 
power  and  configuration. 

On  the  other  hand,  global  efficiency  is  defined  as  the  ratio  of  the  real 
control  power  of  a  globally  optimal  control  based  on  a  spatially  continuous  input 
profile  to  the  real  control  power  of  achieving  the  same  control  objectives  with 
a  spatially  discontinuous  input  profile.  The  definition  is  based  on  dynamically 
similar  control  systems  (Ref.  1).  Physically,  the  global  efficiency  is  related 
to  the  degree  of  extra  control  power  associated  with  a  spatially  discontinuous 
input  configuration  in  controlling  a  DPS  versus  accomplishing  the  same  task  with 
a  spatially  continuous  input  profile  and  therefore,  it  should  be  a  predominant 
indicator  of  the  effectiveness  of  the  input  configuration  to  control  DPS. 

Modal  efficiency  coefficient  is  a  proportionality  constant  between  the 
global  efficiency  and  the  relative  model  efficiency.  Thus,  with  these  quan¬ 
tities,  the  effect  of  finite  dimensionalization  of  the  control  design  model  and 
the  effect  of  nature  of  spatial  discretization  of  the  input  profile  in  con¬ 
trolling  a  DPS  can  be  studied. 

The  concept  of  efficiency  is  introduced  to  address  particular  issues  of  DPS 
control.  A  good  physical  understanding  of  the  concept  is  best  brought  about  by 
considering  a  DPS  formulation.  Therefore,  the  basic  aspects  of  the  concept  will 
be  presented  from  the  perspective  of  a  DPS  formulation.  On  the  other  hand,  in 
practice,  for  a  complex  structure  only  a  spatially  discrete  formulation  is 
explicitly  available.  Since  such  a  spatially  discrete  formulation  acts  as  a 
surrogate  for  DPS  formulation,  the  transition  of  the  definitions  of  efficiency 
to  spatially-discrete  structural  dynamics  is  also  presented.  In  particular, 
since  spatially  discrete  models  are  almost  always  obtained  via  the  Finite  Element 
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method  (FEM) ,  specific  interpretations  of  the  concepts  for  the  finite  element 
models  are  given. 

The  premise  of  the  concepts  presented  in  this  report  is  that  a  structure- 
control  system  must  strive  to  achieve  the  highest  possible  efficiencies.  In 
that,  the  control  system  should  be  designed  such  that  as  much  of  the  real  control 
power  as  possible  should  be  channeled  to  the  reduced-order  control  design  model 
leaving  little,  or,  of  possible,  no  control  power  spillover  for  the  truncated 
dynamics.  This  objective  requires  maximization  of  the  relative  model  efficiency. 
Furthermore,  any  control  design,  characterized  by  a  spatially  discontinuous  input 
profile,  should  try  to  approximate  the  performance  of  the  globally  optimal  solu¬ 
tion  characterized  by  a  dynamically  similar  spatially  continuous  input  profile. 
This  object  requires  maximization  of  the  global  efficiency.  Based  on  the  pro¬ 
posed  concepts,  an  efficient  model  reduction  approach  can  be  proposed  in  that  for 
a  given  input  configuration,  the  components  which  contribute  least  to  the 
efficiencies  can  be  discarded. 

The  efficiencies  will  be  shown  to  be  functions  of  control  parameters  such 
as  the  number  of  control  inputs  and  locations,  method  of  control,  and  model 
order;  as  well  as  structural  parameters  such  as  natural  frequencies  and  mode 
shapes.  All  relevant  parameters  can  be  altered  to  change  the  efficiency  of 
hybrid  structure-control  system.  An  important  aspect  of  the  results  is  that  the 
performance  of  the  entire  infinite-dimensional  (®-D)  DPS  can  be  studied  based  on 
quantities  computed  from  the  finite-dimensional  control  design  model  alone 
without  any  knowledge  of  the  truncated  dynamics.  The  concepts  introduced  are 
valid  irrespective  of  the  theory  or  method  by  which  the  controls  are  designed. 
The  efficiency  concept  can  be  used  both  for  analysis  and  design  of  the  structure - 
control  system. 
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2.0  EFFICIENCY  OF  STRUCTURE -CONTROL  SYSTEMS 


2 . 1  Introduction 

The  concept  of  efficiency  for  structure-control  systems  is  introduced  in 
this  chapter.  The  developments  presented  in  this  chapter  constitute  the  central 
foundation  and  concepts  around  which  this  investigation  revolves.  A  model 
efficiency  and  a  global  efficiency  are  defined  for  the  structure -control  system 
from  the  point  of  view  of  control  power  consumed  during  control.  The  efficiency 
concept  is  illustrated  to  be  a  useful  tool  in  understanding  the  interaction 
between  the  control  variables  and  structure  variables  and  is  shown  to  provide 
insight  to  the  behavior  of  the  structure -control  system.  This  chapter  and  the 
illustrations  of  the  theory  presented  herein  focus  on  the  analysis  of  a 
structure -control  system.  The  illustrative  examples  evaluate  performances  of 
several  linear  quadratic  regulator  designs  for  the  ACOSS-4  structure  (Ref.  10) . 

In  Section  2.2,  control  objectives  for  a  distributed  parameter  system  (DPS) 
are  stated.  Section  2.3  defines  a  global  control  power  for  the  DPS  and  Section 
2.4  presents  the  globally  optimal  control  for  the  DPS.  Suboptimal  control  of  the 
DPS  is  discussed  in  Section  2.5.  Section  2.6  generalizes  the  concepts  of  control 
power  to  spatially  discrete  systems,  in  particular  finite  element  models  (FEM) 
of  structural  systems  are  noted.  The  definition  of  efficiency  of  a  structure- 
control  system  is  introduced  in  Section  2.7.  Also,  presented  in  Section  2.7  is 
the  methodology  for  the  use  of  efficiency  as  an  analysis  and  design  tool. 
Finally  illustrative  examples  are  given  in  Section  2.8  for  the  control  of  ACOSS-4 
tetrahedral  truss  structure. 

2 . 2  Control  Objectives  for  the  DPS 

We  consider  the  DPS  equation  of  motion 

m(p)u(p,t)  ♦  £[u(p, t) ]  =f(p,t)  (2.1) 
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where  in,  u,  £,  and  f  are  mass  distribution,  displacement  field  vector,  stiffness 
operator  matrix,  and  the  input  field  vector,  respectively,  p  denotes  a  position 
vector,  which  will  be  suppressed  in  the  following  for  convenience.  In  general, 
(2.1)  will  represent  the  three-dimensional  partial  differential  equations  of 
motion.  For  convenience,  we  shall  assume  £  >  0. 

The  eigenvalue  problem  of  (2.1)  is 

=  t  4>t  (2.2) 

with  the  orthogonality  relations 

£(p)tfin*.dZ>(p)  -  ,  Jtorf.dD(p)  «  (#„  ,  r ,  s  =  1,2,...  (2.3) 

where  wr  is  the  natural  frequency,  4>r  is  the  corresponding  vector  of  eigen¬ 
functions,  and  D(p)  is  the  structural  domain.  Introducing  the  modal  expansion 

u=E^(P)Ct(C)  (2.4) 

r 

(2.1)  is  transformed  to 

t(c)  +  u?ft(t)  *  f'(t)  r  =  1,2,...  (2.5) 

in  the  modal  configuration  space  where  £r  is  a  modal  configuration  coordinate  and 
fr(t)  represents  a  modal  input  coordinate.  Similar  to  modal  expansion  (2.4)  of 
the  displacement  field  vector,  a  modal  expansion  for  the  input  field  vector  can 
be  written 

f(P>  O  =  Em<P)^i(P)fr(t)  (2.6) 
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From  (3)  and  (6)  the  modal  Input  coordinate  fr  can  be  obtained  as 

fx(t)  =  J^(p)  f(p,t)dD  (2.7) 

The  control  objective  on  the  «-D  DPS  is  to  insure  proper  allocation  of  n 
pairs  of  eigenvalues  of  a  set  of  {n}  modes  from  Equation  (2.5).  Another 
objective  is  to  minimize  control  spillover  effects  so  that  the  response  of  the 
control  design  will  not  be  degraded  by  excessive  control  spillover.  In  improving 
stability  characteristics  of  time- invariant  linear  systems,  the  premise  of  all 
control  methods  is  to  obtain  desirable  eigenvalue  locations  either  directly  or 
indirectly.  The  phrase  "eigenvalue  allocation"  is  used  here  in  a  general  sense 
to  address  both  direct  and  indirect  means  of  allocation.  As  examples,  a  direct 
eigenvalue  allocation  technique  is  the  Simon-Mitter  Algorithm  or  any  other  tech¬ 
nique  in  which  desired  eigenvalue  locations  are  imposed  as  explicit  constraints 
in  obtaining  the  control  gains;  on  the  other  hand,  an  indirect  eigenvalue  alloca¬ 
tion  technique  is  to  use  the  linear  quadratic  regulator  (LQR)  theory  by  which 
eigenvalue  positions  are  obtained  indirectly  as  a  by-product  of  optimizing  a 
performance  measure.  Hence,  the  LQR  approach  is  sometimes  categorized  as  an 
optimal  eigenvalue  allocation  technique.  No  implication  is  intended  in  this 
paper  as  to  the  use  or  necessity  of  direct  eigenvalue  allocation  techniques  in 
understanding  and  applying  the  proposed  concepts  of  efficiency.  In  the  sequel, 
the  emphasis  will  be  on  the  generic  meaning  of  the  phrase  "eigenvalue  allocation" 
as  a  qualifying  phrase  for  the  function  of  the  control  system  which  may  have  been 
designed  by  either  a  direct  or  indirect  allocation  whatever  the  case  may  be. 

From  the  control  objectives  stated  above,  one  would  infer  the  ideal  control 
to  be  the  one  by  which  (n)  pairs  of  closed-loop  eigenvalues  are  located  as 
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desired  with  the  minimum  control  power,  and  control  spillover  is  eliminated  com¬ 
pletely.  The  solution  for  such  an  ideal  control  for  the  ®-D  DPS  (2.1)  will  be 
stated  in  Sec.  (2.4). 

2 . 3  Global  Control  Power  for  the  DPS 

Whatever  one's  favorite  control  design  may  be  to  satisfy  the  control  objec¬ 
tives,  one  can  define  and  compute  a  global  control  evaluation  functional  for  the 
°°-D  DPS  in  the  form 

<o 

S  =  J  JVl(p)fT(p,t)f(p,t)dD  dt  (2.8) 


(2.8)  represents  the  total  quadratic  control  power  expended  on  the  actual  DPS. 
The  control  power  (2.8)  is  dimensionally  the  average  power  consumed  by  the 
control  design  over  the  control  period.  Because  (2.8)  is  the  control  power  on 
the  entire  DPS,  it  is  recognized  as  a  global  quantity.  Since  the  global  control 
power  defined  by  (2.8)  is  computed  from  the  physical  input  field  vector,  it  will 
also  be  referred  to  as  the  real  control  power  SR  where  superscript  R  denotes 
"real,"  hence 

SR  =  JjV^pK^p.t)  f  (p,t)  dD  dt  (2.9) 

SR  can  be  computed  for  any  given  f(p,t)  regardless  of  the  details  of  the  control 
design  technique  by  which  it  is  computed.  By  a  similar  motivation,  we  define  a 
global  modal  control  power  functional 
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(2.10) 


f  fr2(t)dt 

Jo 

where  M  denotes  that  the  quadratic  controls  are  the  modal  control  input  coordi¬ 
nates  and  the  subscript  ®  implies  that  all  modes,  hence  the  entire  DPS,  are  con¬ 
sidered  in  the  computation.  In  contrast  to  SR,  which  involves  the  real  input 
field,  SM  seems  to  represent  an  abstract  quantity  since  modal  inputs  are  used  in 
its  definition.  However,  SR  and  S3  are  related. 

Substituting  the  modal  expansion  (2.7)  into  definition  (2.9)  and  using  the 
orthogonality  relations  one  obtains 

SR  =  m-Ma  dD  dt  -  Y.  f^r^)  dc  (2.11) 

J J  r « 1  8 ■ 1  r*l J 

hence 

SR  =  Si  (2.12) 

Identity  (2.12)  defines  an  invariance  property.  The  global  control  power 
for  the  DPS  is  frame  indifferent;  it  is  the  same  whether  one  studies  it  in  the 
real-space  of  (2.9)  or  the  modal-space  of  (2.10). 

The  global  modal  control  power  S3  can  be  decomposed  into 

n  •  m 

sS  =  £  f  fr2(t)  dc  ♦  Y  fof*2(t)  dt  (2.13) 

r*l J  0  r *n+l J 


r-1 
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s%  =  s*  +  s£ 


(2.14) 


and  from  (2.12) 

SR  =  S*  +  S%  (2.15) 

where  the  definitions  of  and  Sjj  should  be  evident  from  (2.13)  and  (2.14).  S” 
is  the  portion  of  the  control  power  SR  channeled  for  control  of  the  {n}  modes  and 
S”  is  the  remaining  control  power  channeled  into,  or  better  said,  spilled  over 
to  the  uncontrolled  modes.  If  the  set  of  modes  (n)  are  referred  to  as  the  con¬ 
trol  design  model,  s£  will  be  termed  design  model  control  power.  Similarly,  S™ 
will  be  referred  to  as  control  power  spillover. 

Since  SR,  SfJ,  S”  are  positive  definite  quantities: 

SR>S“  (2.16) 

where  the  equality  is  satisfied  if  and  only  if  S^J  -  0;  that  is,  when  there  is  no 
control  power  spillover. 

An  important  feature  of  the  control  of  the  <»-D  DPS  is  imbedded  in 
inequality  (2.16).  Because  SR  is  computed  by  using  the  real  input  f(p,t)  applied 
to  the  actual  structure,  and  S”  is  computed  by  using  the  modal  inputs  to  the 
finite-dimensional  control  design  model,  inequality  (2.16)  relates  how  the 
control  design  model  performance  stands  relative  to  the  actual  DPS.  It  is  clear 
that  any  mismatch  between  SR  and  S^  would  automatically  mean  that  some  of  the 
control  power  is  lost  to  the  residual  modes,  Sy  *  0.  On  the  other  hand,  in 
accordance  with  the  control  objectives,  the  ideal  control  system  for  the  DPS 
would  yield  Sjj  —  0 ,  that  is,  it  would  minimize  the  power  spillover 
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It  remains  to  address  the  specifics  of  how  one  might  realize  a  minimum 
global  control  power  to  achieve  the  control  objectives.  To  this  end,  we  shall 
assume  linear  state- feedback  control. 

2 . 4  Globally  Optimal  Control  for  the  DPS 

In  the  absence  of  other  objectives  and  design  constraints,  it  is  reasonable 
to  try  to  achieve  the  control  objectives  stated  in  Sec.  (2.2)  with  the  minimum 
amount  of  global  control  power.  Hence  one  can  state  the  optimization  problem 
Minimize  S5,  or  equivalently,  SR 
subject  to: 

p{±  iut)  -*  p{at  ±  if 3t)  r  =  1,2 . ®  (2.17) 

This  is  an  optimization  problem  for  the  ®-D  DPS.  The  uncontrolled  eigen¬ 
values  ±iwr  are  relocated  to  specified  positions  ax±ipr  where  />!)  represents  an 
eigenvalue  spectrum.  There  are  no  restrictions  on  the  constraint  values  aI  and 
f)z.  The  case  where  only  n  pairs  of  eigenvalues  are  relocated  to  new  positions 
and  the  remaining  residual  pairs  are  not  moved,  is  a  special  case  of  the  above 
formulation  since  we  can  always  write 

p{±iwr}  -»  p{a,  ±  r  =  1,2, .  . .  ,n 


pl+iWj)  -*  p{±iwt)  r  *  n+1 . « 


(2.18) 


The  solution  of  the  minimum  global  control  power  is  known  to  be  (Refs.  1,2) 
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-  (&v;(o  -  &*2c  ;(t)} 


(2.19) 


r 


1.2 _ 


aj 


St 


=  W, 


- 


Si. 


2ar 


(2.20) 


where  (*)  denotes  the  optimum  quantities.  Substituting  the  solution  (2.19)  into 
(2.6)  one  obtains  the  optimum  input  field  vector  which  can  achieve  the  desired 
eigenvalue  locations  with  minimum  control  power  (Ref.  3): 

f(p,t)  -  J((?i(p,p')  u(p',t )  +Gl(p,p')  u(p',t)dD(pf)  (2.21) 

where  G*  and  G2  are  identified  as  symmetric  optimum  distributed  control  influence 
(Kernel)  functions 

CO 

GI(P.P')  =  Y,  Sr*inj(p)^r(p)^r(P/)iD(P/)  =Gi‘(P/.P)  (2.22) 


Gz(p,p')  =  X]  St2nj(P)<Ji,r(P)^r(p/)m(p/)  =  Gl(p'  ,p) 


We  observe  the  following  characteristics  of  the  optimal  control  solution 
(2.21):  the  optimal  modal  input  coordinate  fr*  is  a  feedback  of  only  the  corre¬ 
sponding  r-th  modal  coordinate,  therefore  optimal  modal  control  coordinates  are 
independent  of  each  other.  This  feature  of  feedback  control  has  come  to  be  known 
as  independent  modal-space  control  (Ref.  4).  The  corresponding  optimal  input 
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field  f*(p,t)  is  spatially  continuous  since  the  modal  synthesis  of  spatially  con¬ 
tinuous  functions  m(p)^r(p)  is  a  spatially  continuous  function.  Without  proof, 
we  also  state  that  controlled  DPS  under  the  optimal  spatially  continuous  feedback 
input  has  the  same  eigenfunctions  as  the  uncontrolled  DPS,  preserving  its  natural 
properties  (Refs.  1,2).  Therefore,  the  optimal  control  (2.19-2.21)  has  also  been 
referred  to  as  Natural  Control  (Ref.  2). 

Specifically,  if  the  desired  closed-loop  eigenvalues  are  given  as  the  set 
(2.18)  from  (2.19,  2.20)  we  compute  g*x  -  g*2  -  0  for  r-n+1,  ...,  «>  which  yields 
f  j  —  0  ,  r-*n+l ,  n+2 

Upon  substituting  this  result  into  (2.21)  we  get 

f*(p,t)  =  J{G1*n(p,p/)u(p/,t)  +  (?2n(P>P/ )  u(p',t)}  dD(p')  (2.23) 

where  Gjn  and  G2n  are  the  same  as  (2.22)  except  that  the  summation  ends  at  n.  The 
point  is  that  the  solution  (2.23)  does  not  represent  a  model  truncation,  instead 
the  required  summations  end  at  n  because  the  remaining  terms  have  been  computed 
to  be  zero. 

The  above  procedure  indicated  that  optimal  control  can  be  found  with  vir¬ 
tually  no  effort  for  any  eigenvalue  set  for  the  ®-D  DPS.  It  remains  to  check  the 

control  spillover  effect  over  the  uncontrolled  modes  r-n+1 .  To  this  end, 

we  substitute  the  form  of  f*(p,t)  given  by  (2.23)  into  (2.7)  which  yields  upon 
recognizing  the  first  orthogonality  in  (2.3), 

fs(t)  =  0  s  =  n+1 ,  .  . 

Perfect  spillover  elimination  from  residual  modes  is  also  achieved  by  the  optimal 
control  (2.23).  We  must  point  out  that  the  spillover  inputs  fs(t)  vanish  due  to 
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appearance  of  <£r(p)  in  f(p,t)  regardless  of  the  functional  form  of  the  modal  con¬ 
trol  inputs  fr(t),  r-1,2 . n;  that  is  whether  modal  inputs  are  independent  or 

not.  Therefore,  spillover  control  is  ultimately  not  a  matter  of  what  the  tem¬ 
poral  behavior  of  the  control  inputs  is,  but  is  a  matter  of  spatial  distribution 
of  control.  Any  other  spatial  distribution  of  input  would  not  yield  perfect 
spillover  elimination  at  least  theoretically. 

Last  but  not  the  least  important  feature  of  the  optimal  control  is  that  the 
solution  is  unique  (Refs.  1,2),  therefore  it  is  globally  optimal,  and  it  controls 
the  ®-D  DPS  accomplishing  the  control  objectives  ideally. 

The  global  control  power  for  the  optimal  control  for  the  DPS  can  be  evalu¬ 
ated  by  substituting  (2.21)  or  (2.23)  into  (2.8) 

S*  =  J/n-1f*T(p,t)f*(p,t)dD  (2.24) 


which  yields  from  (2.12) 


S 


R  _ 


=  s* 


(2.25) 


2 . 5  Subontimal  Control  for  the  DPS 

By  definition,  any  control  input  field  f(p,t)  of  the  form 
f (p, t)  *  f*(p, t) 

will  have  a  higher  real  control  power  than  S*.  The  most  common  suboptimal 
control  is  the  one  that  seems  most  practical  to  implement;  it  is  the  point  (or 
localized)  input  distribution 

m 

=  E  *(p-APk)Fk(t)  (2.26) 

k»  1 
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where  S  is  the  spatial  Dirac  delta  function  Apk  is  the  domain  of  the  influence 
of  the  local  k-th  input  and  Fk(t)  is  the  total  input  over  Apk. 

The  real  control  power  for  the  suboptimal  control  is  again  computed  by 
using  (2.26)  in  (2.9).  Denoting  the  total  input  vector  of  m  inputs  by  F=[Fi 
F2...Fm]T,  it  can  be  shown  that  (Ref.  5) 

5r  =  J>(t)TR  F(t)  dt  i?  =  diag[m(pk)Apk]'1  (2.27) 

and  from  identity  (2.12)  the  total  modal  control  power  corresponding  to 
suboptimal  control  profile  (2.26)  is 

S“  =  j>T.R  F  dt 

Furthermore,  the  design  model  control  power  for  {n}  pairs  of  relocated 
eigenvalues  is 

Sc  =  £  k2(t)  dt  -  ^  ( t )  f n  ( t)  dt  =  fFTflX^  dt  (2.28) 

I>1 J  J  J 


where  fn(t)  is  the  n- component  design  model  modal  input  vector  generated  by  the 
control  (2.26)  and 

fa(t)*BnF  Bnrk  =  [*rk<ft>]  r  =  1 , 2 . n 

k-1,2,  .  .  .  ,m  ^>rk(pk)  is  the  area  under  <f>k  over  Apk. 

The  control  power  spillover  due  to  localized  inputs  F  can  be  evaluated  as 

CO 

s,7  =  Jfr2(t)dt  »  JfJ(t)fn(t)dt  =  Jf^Xf  dt  (2.29) 
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where 


fu(t)  =  BirF ,  Byik  =  [ <^>rk (Pk)  ]  r  =  n  +  i . k  =  1 , 2  .  .  .  ,ra 

Bn  and  Ba  are  the  respective  partitions  of  the  B  matrix  for  the  n-th  order  con¬ 

trol  design  model  and  the  uncontrolled  dynamics. 

One  does  not  need  to  compute  the  control  spillover  performance  according 
to  (2.29)  if  SR  and  S”  are  already  available.  For,  from  (2.27)  and  (2.28) 

S“  =  SK  -  S“  =  -  B*Bn]  F  dt  (2.30) 

in  which  F(t)  is  directly  available  from  the  control  design  model  once  it  is 

selected.  Equation  (2.30)  indicates  that  for  a  given  DPS  (hence  the  mass  dis¬ 

tribution  and  therefore  the  matrix  R  are  known)  when  a  control  design  model  is 
selected,  the  control  spillover  performance  can  be  determined  solely  on  the  basis 
of  the  control  design  model.  No  knowledge  of  the  uncontrolled  dynamics  is 
needed.  This  points  out  the  usefulness  of  the  judicious  definition  of  the  con¬ 
trol  evaluation  functional  S  in  the  form  of  (2.9). 

The  control  power  spillover  for  the  suboptimal  point  input  profile  (2.26) 
cannot  vanish,  therefore  from  (2.15)  we  deduce 
SR  >  s£ 

and  since  SR  as  given  by  (2.27)  is  suboptimal 
SR  >  S* 

2 . 6  Control  Powers  for  D i sc  re te - Svs terns 

Quite  often,  instead  of  partial  differential  equations,  a  large  set  of 
spatially  discretized  ordinary  differential  equations  are  assumed  to  describe  the 
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dynamics  of  the  DPS  adequately,  such  as  the  finite  element  models  (FEM)  of  com¬ 
plex  structural  systems.  Extension  of  the  previous  definitions  and  results  to 
such  cases  will  be  useful.  We  assume  that  the  structural  system  is  described  bv 

Mq(£)  +  Kq(t)  =  Q(t)  (2.31) 

instead  of  (2.1),  where  q(t)  and  Q(t)  are  each  N-component  generalized  coordi¬ 
nates  and  the  generalized  forces  vectors,  respectively.  M  and  K  are  N  x  N  sym¬ 
metric,  positive  definite  mass  and  stiffness  matrices,  respectively.  Equation 
(2.31)  is  usually  referred  to  as  the  N-dimensional  evaluation  model  replacing  the 
®-D  DPS,  (2.1)  . 

Denoting  by  E  the  modal  matrix  associated  with  the  system  (2.31),  the  modal 
transformations  and  the  orthogonality  relations 

g  =  ££  ,  f ( t )  *  ErQ ( t )  ETME  =  I  ,  EtKE  =  [w2]n  (2.32) 

again  yield  the  modal  equations  of  motion  (2.5)  with  the  exception  that  this  time 

r=l,2 .  N  represents  the  complete  system.  Again,  by  definition  the  total 

modal  control  power  is 

3?  ■  E  dt  t)  f(t)  dt  (2.33) 


where  N  denotes  the  total  system  as  opposed  to  <*>  in  (2.13)  for  the  DPS  Equation 
(2.1).  The  corresponding  real  control  power,  after  recognizing  the  invariance 
property  (2.12)  is 

SR  =  Sg  =  f fTf  dt  *  [QTEETQ  dt  (2.34) 
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Noting  from  the  orthogonality  relations  that  EET  -  M'1,  for  a  discrete-system  the 
equivalent  of  SR  in  (2.9)  is 

SR  =  JqtM_1Q  dt  (2.35) 

Hence,  the  generalized  input  vector  Q(t)  plays  the  role  of  f(p,t). 

In  particular,  if  (2.31)  represents  the  FEM  equations  of  motion,  the  gen¬ 
eralized  loads  vector  Q(t)  is  the  vector  of  joint  loadings.  If  F(t)-[F1  F2  .  .  . 
Fn]t  is  the  real  joint  inputs  vector  of  f(p,t)  one  can  write 

Q  =  DF  (2.36) 

where  D  is  the  joint  loads  distribution  matrix.  If  f(p,t)  is  spatially  con¬ 
tinuous  over  the  whole  structural  domain,  then  by  necessity,  it  will  yield  an 
input  at  each  joint  along  every  joint  degree  of  freedom.  It  follows  that  the 
equivalent  of  a  spatially  continuous  input  f(p,t)  in  a  FEM  setting  is  tantamount 
to  having  a  full  generalized  loads  vector  Q.  On  the  other  hand,  if  the  input 
f(p,t)  is  spatially  discontinuous,  such  as  in  (2.26),  that  will  be  tantamount  to 
having  a  Q  with  some  zero  components.  If  there  are  m  independent  joint  inputs 
F1(  ...  Fm  as  in  (2.26)  substitution  of  (2.36)  into  (2.34)  yields 

SR  =  ^F(t)TR  F(t)  dt  R  =  DrM~1D  (2.37) 

Hence,  the  weighting  matrix  R  in  (2.37)  is  the  FEM  equivalent  of  the  weighting 
matrix  R  in  (2.27)  for  the  partial  differential  equations  of  motion. 

It  is  clear  from  the  form  of  R  that,  given  any  F(t)  ,  SR  describes  the 
control  power  for  the  entire  evaluation  model.  Even  if  F(t)  may  have  been 
designed  by  considering  only  a  reduced  {n)th  order  modal  model  of  the  N-th  order 
system  (2.31),  when  SR  is  computed,  it  will  be,  according  to  (2.37),  the  global 
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control  power  for  the  total  N-th  order  evaluation  model,  not  for  the  n-th  order 
reduced-control  design  model. 

For  an  n-th  order  control  design  model,  it  is  easy  to  see  that  the  counter¬ 
parts  of  (2.13)  and  (2 . 27) - (2 . 29)  are 

S“  =  S“  +  s£  (2.38) 

where 

S *  =  j>TDT£n££  D  F  dt  ,  S%  =  jVTDTEu£^  D  F  dt  (2.39) 

S%  •  SR  -  S“  =  [FTDT[/rx  -  FnFj]D  F  dt  (2.40) 

and  E„  and  Eu  are  the  control  design  model  and  the  uncontrolled  model  modal 
matrices,  respectively.  Here  again,  from  (2.40),  for  a  given  physical  system 
(hence  the  evaluation  model  mass  matrix  M  is  known)  and  the  n-th  order  control 
design  model,  control  spillover  performance  can  be  ascertained  solely  on  the 
basis  of  the  control  design  model.  Specifically,  if  a  FEM  is  used,  the  modes 
that  are  poorly  computed  will  be  inconsequential  from  the  control  point  of  view 
as  long  as  those  modes  are  ir  the  uncontrolled  set. 

Finally,  it  remains  to  ascertain  the  counterparts  of  the  globally  optimal 
control  and  the  control  power  (2 . 23) - (2 . 25)  for  the  system  of  (2.31).  From  the 
previous  discussions  and  the  nature  of  the  globally  optimal  control,  it  is 
straight  forward  to  obtain 

S*  =  f(?*TW1Q*dt  »  fF,TD*Ttf-1D*F*dt  (2.41) 
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Q*  =  M  E  f*(t)  =  D’F*  ,  F*  =  D-'M  E  f*(t) 


(2.42) 


where  f*(t)  is  the  vector  of  N- independent  modal  inputs  fr*  r-1 , 2 ,  .  .  .  ,N  precisely 
as  computed  according  to  Equations  (2.19),  (2.20).  For  a  FEM,  Q*  is  the  full 
generalized  inputs  vector  tantamount  to  having  N  joint  inputs  Fr*,  r-l,2,...,N, 
that  is  as  many  inputs  as  the  total  number  of  degrees  of  freedom.  In  other 
words,  D*  for  the  globally  optimal  control  must  be  a  full  rank  N  joint  loads 
distribution  matrix. 

The  forms  of  Q*  and  its  interpretation  makes  it  clear  that  a  FEM  not  having 
inputs  at  each  joint  along  each  joint  degree  of  freedom  would  correspond  to  a 
discontinuous  input  profile  and  have  a  suboptimal  performance. 

2 . 7  Efficiency  of  a  Structure-Control  System 

Implementable  control  designs  for  large  flexible  structures  such  as  complex 
truss- like  configurations  planned  for  the  space  station  will  inevitably  employ 
spatially  discontinuous  suboptimal  input  profiles  consisting  of  a  large  number 
of  distributed  point  force  and  torque  actuators.  In  view  of  the  control  objec¬ 
tives  and  the  features  of  the  control  powers  we  discussed  heretofore,  it  would 
be  desirable  for  any  implementable  structural  control  system  to  channel  as  much 
of  the  real  control  power  as  possible  to  the  control  design  model.  In  other 
words,  the  power  spillover  S”  should  be  minimized  by  the  control  design.  An 
equally  desirable  feature  of  the  control  design  would  be  to  keep  the  total 
control  power  as  small  as  possible,  that  is,  to  keep  the  real  control  power  of 
the  design  as  near  to  the  globally  optimal  control  power  as  possible.  These 
aspects,  by  necessity,  bring  about  the  concept  of  efficient  structure-control 
designs.  An  important  element  of  the  structure-control  design  process  must  be 
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to  find  the  most  efficient  structure -control  combination  for  the  control 


objectives . 

We  define  the  percent  global  efficiency  of  a  structure -control  system  as 

e*%  =  —  x  100  <  100%  (2 .43) 

SR 

where  cR  is  the  real  control  power  of  any  suboptimal  control  design  with  the 
clostd-loop  eigenvalue  spectrum  {p)-{p).  S*  is  the  globally  optimal  control 
power  corresponding  to  the  same  eigenvalues  (cost  0).  This  is  to  say  that  the 
global  efficiency  is  based  on  the  comparison  of  dynamically  similar  (Ref.  1) 
globally  optimal  and  suboptimal  control  designs  for  the  desired  closed- loop 
eigenvalues.  Since  SR  is  suboptimal,  the  upper  bound  of  global  efficiency  is 
100%. 

Next,  we  define  the  percent  relative  model  efficiency  of  a  structure - 
control  system  to  be 

SM 

e%  =  _£  x  100  <  100%  (2.44) 

SR 


The  relative  model  efficiency  is  an  indicator  of  the  percentage  of  the  real 
control  power  channeled  into  the  control  design  model,  the  balance  indicating  the 
control  power  spillover.  This  efficiency  is  determined  solely  by  using  the 
properties  of  the  particular  control  design  model.  There  is  no  reference  to  the 
corresponding  globally  minimum  control.  Therefore,  we  refer  to  e  as  the  relative 
model  efficiency.  A  less  than  perfect  model  efficiency  automatically  implies 
control  power  wasted  to  uncontrolled  modes.  However,  as  a  100%  e  means  no 
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control  spillover,  it  will  not  guarantee  a  100%  global  efficiency  as  SR  and  S* 
may  still  be  different. 

The  relative  model  efficiency  and  the  global  efficiency  are  related  by 


e  =  ne* 


(2.45) 


where  n  is  defined  to  be  the  modal  efficiency  coefficient. 

Complimentary  to  the  above  definitions,  one  can  also  introduce  the  global 
and  the  model  spillover  quotients 


sq* 


$ 

S* 


(SR-S“)  _  1 

S*  e* 


=  1  -  e  <  1 


(2.46) 


sq  indicates  the  portion  of  the  real  control  power  lost  as  control  spillover 
cost.  sq*  indicates  the  control  power  spillover  of  the  suboptimal  design  as  a 
fraction  of  the  globally  minimum  control  power  that  would  be  expended  on  the 
entire  DPS.  Studies  show  that  the  control  power  spillover  a  suboptimal  control 
profile  can  incur,  can  be  many  times  more  than  it  would  take  to  control  the 
entire  system  with  a  spatially  continuous  optimal  input. 

Given  an  initial  modal  state  disturbance  xQ  for  a  stable  control  system, 

x0  =  [^(0)  ^(0)  .  .  .  £n(0Kn(0)] 


The  control  power  for  infinite  time  control  are  given  by 


SR  =  XqPrx0  ,  S“  =  XqP*xq  ,  S *  =  XqP*x0 


(2.47) 
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where  PR,  P*1  and  P*  are  the  real,  modal,  and  the  globally  optimal  (natural)  con¬ 
trol  power  matrices.  PR  and  PM  can  be  obtained  as  the  solutions  of  the  asso¬ 
ciated  Lyapunov  equations  for  any  suboptimal  control  design  discussed  in  the 
preceding  sections.  The  relevant  equatioins  are  presented  in  Chapter  3.  The 
natural  control  power  matrix  P*  is  given  in  closed  form  in  Ref.  (6)  and  in 
Appendix  B. 

The  global  and  relative  model  efficiencies  of  any  control  design  can  now 
be  computed  by  using  the  cost  matrices 


XqP*x0  XqP“x0  XqP”xq 

XqPrx0  XqPrx0  XqP’Xq 


(2.48) 


Each  one  of  the  efficiencies,  through  the  power  matrices,  depends  on: 

*  The  number,  type  and  locations  of  localized  inputs  Fk(t) 

*  The  particular  control  design  technique  used  to  compute  the  actual 
spatially  discontinuous  feedback  input  F 

*  The  order  n  of  the  control  design  model  and  the  closed- loop  eigen¬ 
value  spectrum  (p) 

*  Structural  parameters  through  the  appearance  of  modal  frequencies  and 
mode  shapes 

*  The  initial  modal  disturbance  state  x0. 

For  the  analysis  and  design  of  structural -control  systems  via  efficiencies, 
one  would  typically  take  the  following  steps:  1)  For  any  set  of  selected  system 
variables  and  parameters  mentioned  above  (such  as  a  given  n-th  order  design  model 
and  a  control  input  configuration)  obtain  a  control  law  by  whatever  technique  or 
theory  deemed  appropriate.  2)  For  the  control  inputs  obtained  in  step  1,  compute 
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the  corresponding  SR  and  s£  defined  by  Eqs .  (2.27),  (2.28)  or  Eqs.  (2.37)  and 
(2.39),  and  calculate  the  relative  model  efficiency  e  and  the  model  spillover 
quotient  sq  as  defined  by  Eqs.  (2.44)  and  (2.46),  respectively.  3)  For  different 
values  and/or  sets  of  variables  and  parameters  repeat  steps  1  and  2,  compare  the 
corresponding  model  efficiencies,  simulate  if  necessary,  and  identify  satisfac¬ 
tory  designs  consistent  with  the  designer's  criteria  and  constraints.  In  apply¬ 
ing  steps  1-3,  one  should  recognize  that  there  is  no  need  for  explicit  knowledge 
of  the  closed-loop  eigenvalues  if  studies  based  on  relative  model  efficiencies 
are  all  that  is  desired.  However,  in  addition,  if  global  efficiencies  e*  and 
global  spillover  quotients  sq*  are  also  desired  for  further  consideration,  one 
must  then  proceed  with  the  following  steps:  4)  For  the  controls  designed  in  step 
1,  compute  the  corresponding  closed- loop  eigenvalues  { p }  if  they  are  not  alrea  dy 
available.  Otherwise,  this  step  is  not  needed.  5)  For  the  spectrum  (p)  found 
in  step  4,  compute  the  modal  control  gains  and  the  modal  inputs  given  by  Eqs. 
(2.19)  and  (2.20)  of  the  globally  optimal  spatially  continuous  control  which  is 
dynamically  similar  to  the  control  design  of  step  1.  6)  In  accordance  with  Eqs. 
(2.25),  (2.11)  and  (2.12)  obtain  the  globally  minimum  real  control  power  S* 
possible  for  the  eigenvalue  spectrum  {pi  elicited  by  the  control  design  of  step 
1.  Closed-form  solution  for  S*  is  given  in  Appendix  B  for  any  defined  { p ) .  7) 
Calculate  the  global  efficiency  e*  and  global  spillover  quotient  sq*  of  the  con¬ 
trol  design  of  step  1,  as  defined  by  Eqs.  (2.43)  and  (2.46)  by  using  SR  from  step 
2  and  S*  from  step  6.  If  desired,  computed  the  model  efficiency  coefficient 
defined  in  Eq .  (2.45).  8)  For  different  values  and/or  sets  of  variables  and 
parameters  repeat  steps  4-7,  compare  the  corresponding  global  efficiencies  and 
spillover  quotients,  simulate  if  necessary,  and  identify  a  satisfactory  designs. 
9)  Study  the  results  of  steps  3-8  collectively  to  evaluate  the  control  designs. 
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The  efficiency  approach  to  structure  control  will  liberate  the  engineer 
from  the  need  of  detailed  knowledge  of  the  unmodeled  modes.  Because  the  behavior 
of  an  infinite  dimensional  system  can  be  studied  and  understood  by  means  of 
computing  its  efficiencies,  which  require  explicit  knowledge  of  only  the  finite 
number  of  modeled  control  modes.  This  feature  should  make  the  efficiency 
approach  to  control  design  a  practical  tool. 

With  such  an  approach,  it  is  possible  to  determine  the  optimal  control 
input  distribution  and  even  the  optimum  eigenvalue  distribution  for  a  given  nth- 
order  control  design  model.  For  a  given  input  field,  the  efficiencies  can  be 
used  to  determine  the  order  n  of  a  control  design  where  model  orders  that  yield 
high  efficiencies  can  be  selected.  In  addition,  because  the  efficiencies  are 
dependent  on  the  input  distribution,  closed- loop  eigenvalues,  and  other  struc¬ 
tural  and  control  parameters,  different  order  control  models  could  become  more 
efficient  simply  by  changing  the  structure -control  system's  configuration  and 
parameters.  In  all  of  these,  the  objective  then  should  be  to  maximize  the  global 
and  relative  efficiencies. 

In  particular,  for  N-th  order  discrete  evaluation  models,  S„  represents  the 
control  power  consumed  in  controlling  n<N  modes,  while  SR  represents  the  power 
consumed  by  all  N  modes.  Therefore,  the  relative  efficiency  e  becomes  a  valid 
nondimens ional  measure  of  the  effects  of  model  order  reduction.  A  similar  state¬ 
ment  holds  true  for  global  efficiency.  Based  on  these  observations,  a  closed- 
loop  Efficient  Model  Reduction  Technique  can  be  formulated  in  that  one  can 
propose  to  retain  in  the  control  design  model  the  modes  to  which  the  relative 
model  and/or  global  efficiencies  ar  most  sensitive  for  any  given  input 
configuration. 
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In  contrast  to  design,  efficiencies  can  be  used  to  evaluate  the  merits  of 
a  given  control  design  since  they  reflect  the  effects  of  many  variables  of  the 
control  problem.  In  this  chapter,  we  shall  illustrate  the  analyses  of  some 
Linear  Quadratic  Regulator  (LQR)  control  designs  based  on  the  efficiency  concept. 
Also,  the  proposed  concept  of  Efficient  Model  Reduction  Technique  will  be 
demonstrated. 


2 . 8  Illustrative  Example:  Control  of  ACOSS-4  Tetrahedral  Truss  Structure 
Analysis  via  Efficiencies 

As  a  demonstration  of  the  utility  of  the  efficiency  concept,  the  perform¬ 
ances  of  various  LQR  control  designs  (step  1)  for  the  ACOSS-4  structure  shown  in 
Fig.  2.1  were  evaluated  for  different  order  modal  control  design  models  and  dif¬ 
ferent  number  of  inputs  (steps  3  and  8).  The  inputs  were  located  at  the  pods  of 
the  structure.  A  twelfth  order  (N“12)  evaluation  model  (2.31),  obtained  via  FEM, 
was  considered.  For  a  given  n-th  order  control  design  and  number  of  inputs  1  < 
ra  <  6,  the  control  designs  were  based  on  the  minimization  of  the  LQR  performance 
measure  (step  1) 


J  =  ^  (xX  x  ♦  FrWc  F )  dt 


Wx  =  block-diac jjjij  ®  -2 


wc  =  [1] 


where  the  LQR  control  design  weighting  matrix  Wc  has  no  relationship  to  the 
weighting  matrix  R  uniquely  defined  for  efficiency  purposes.  The  LQR  design 
approach  essentially  is  an  indirect  eigenvalue  allocation.  Instead  of  requiring 
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explicit  eigenvalue  allocation  one  can  implicitly  admit  the  desired  eigenvalues 
(p)  to  be  those  of  the  LQR  solutions  for  specific  choices  of  weighting  parameters 
q  and  r.  For  each  LQR  steady-state  Riccati  equation  solution,  the  closed-loop 
eigenvalues  were  computed  and  assigned  to  be  the  set  {p}  (step  4)  where  upon  the 
corresponding  globally  optimal  control  power  S*  was  computed  (step  6)  for  the  set 
(p).  For  each  LQR  solution,  SR  and  S”  were  computed  by  solving  the  associated 
Lyapunov  equations  for  the  closed-loop  system  (step  2).  In  the  simulations,  the 
2n-th  order  initial  modal  state  x0  was  assumed  to  be  x0  -  [1  Uj  1  1  wn]  and 
the  uncontrolled  modes  were  initially  undisturbed. 

For  control  design  models  of  order  n~2 ,  4,  6,  8,  10,  12  and  the  input  num¬ 
bers  m-1 ,  2,  3,  4,  5,  6  (steps  3  and  8);  e*,  e,  p,  sq*,  and  sq  for  the  chosen  LQR 
weighting  matrices,  Wx  and  Wc ,  were  computed  (steps  2  and  7).  The  results  are 
shown  in  Figs.  2.2  (steps  3  and  8).  The  model  selections  were  made  by  starting 
from  the  lowest  structural  modes  to  the  higher  ones.  For  brevity,  efficiencies 
and  spillover  quotients  are  tabulated  only  for  n~2  and  n=8  (Tables  2.1  and  2.2). 
Efficiencies  of  other  design  models  and  actuator  configurations  can  be  inferred 
from  the  efficiency  curves. 

From  Figs.  2.2a  and  b,  we  observe  the  interactions  among  the  efficiencies, 
the  order  of  the  control  design  model,  and  the  input  configuration.  For  a  given 
number  of  inputs  the  model  efficiency  increases  with  the  order  of  the  control 
design  model.  However,  this  is  not  necessarily  true  for  the  global  efficiency. 
For  one  and  two  inputs  the  global  efficiency  seems  to  increase  with  model  order, 
but  for  three  inputs  increasing  the  model  order  beyond  n=4  decreases  the  global 
efficiency.  We  also  observe  that  for  a  given  control  design  model,  increasing 
the  number  of  inputs  increases  the  global  efficiency,  but  this  is  not  necessarily 
true  for  the  model  efficiency.  Indeed  for  a  sixth  order  design  model  (n-6), 
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three  or  more  inputs  cause  a  decrease  in  the  model  efficiency.  Therefore,  the 
third,  fourth,  fifth,  and  sixth  actuators  are  located  poorly  with  respect  to  the 
truncated  modes  (n  >  7-12)  such  that  more  of  control  energy  is  lost  as  control 
power  spillover  to  cause  a  drop  in  the  model  efficiency.  A  similar  observation 
is  made  for  1  and  2  inputs  in  regards  to  model  efficiency.  Hence,  it  appears 
that  the  first  3  actuators  represent  a  critical  number  of  inputs  for  this  par¬ 
ticular  structure. 

The  curves  of  sq*  and  sq  would  describe  more  vividly  the  effect  of  model 
truncation.  However,  because  sq*  and  sq  are  related  to  e*,  \i  and  e,  for  brevity 
these  curves  are  rot  shown.  Tables  2.1  and  2.2  list  some  values  of  the  spillover 
quotients.  From  the  sq*  values  given  in  Table  2.1  one  reads,  for  example,  that 
for  m*>l  and  n=2  the  amount  of  control  power  lost  to  model  truncation  is  152  times 
the  total  control  power  that  would  be  required  co  control  the  entire  DPS  with  a 
spatially  continuous  optimal  input  profile  (natural  control).  The  control  powers 
SR,  S”  for  the  LQR  designs  and  the  corresponding  dynamically  similar  natural  con¬ 
trol  powers  S*  for  n=2 ,  8  are  shown  in  Figs.  2.3a  and  b.  The  power  plots  show 
that  natural  control  power  S*  can  be  significantly  lower  than  the  real  power  SR 
obtained  by  using  LQR  designed  controls.  The  distances  among  the  power  curves 
are  indicators  of  the  global  and  model  efficiencies  and  the  modal  efficiency 
coefficient . 

The  LQR  values  of  S8  decrease  monotonically  to  limit  values  with  increasing 
number  of  inputs.  The  dynamically  similar  natural  control  powers  S*  increase 
monotonically  to  different  limit  values.  We  conjecture  that  the  two  designs  will 
not  converge  because  of  fundamental  differences  in  their  design  concepts.  The 
natural  control  is  a  distributed  partial  differential  equation  control  solution 
according  to  Equations  (2.1),  (2.21)  and  (2.22).  On  the  other  hand,  the  LQR 
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solution  is  a  discrete  control  solution  based  on  the  a  priori  reduced-order 
(truncated)  model  of  the  dynamic  system.  In  order  for  the  two  optimal  solutions 
to  con.erge  to  different  limits,  their  conceptual  framework  must  be  inherently 
different.  This  suggests  that  other  than  LQR  discrete  closed- loop  laws  can  be 
formulated  as  direct  approximations  to  the  closed- loop  distributed  natural  con¬ 
trol  solution  with  control  powers  between  the  LQR  and  natural  control  powers 
(Ref.  3). 

The  response  profiles  for  a  sensor  colocated  with  the  first  input  are  given 
in  Fig.  2.4  for  different  design  model  orders  and  inputs.  The  response  profiles 
corresponding  to  the  natural  (globally  optimal)  control  with  continuously  dis¬ 
tributed  input,  suboptimal  control  of  the  n-th  order  control  design  model  with 
m  point  inputs  and  the  evaluation  model,  which  includes  the  control  spillover 
effects  of  the  suboptimal  control,  are  superposed  in  Figs.  2.4a,  b  and  c  for 
comparison  purposes.  It  is  seen  that  almost  identical  responses  can  be  obtained 
with  drastically  different  control  powers.  The  similarity  between  the  responses 
of  the  (suboptimal)  control  design  model  and  the  evaluation  model  for  n-8 ,  ra-2 
in  Fig  2.4b  may  suggest  that  model  truncation  is  insignificant.  This  is  true 
from  an  output  viewpoint.  However,  there  still  exists  a  considerable  ineffi¬ 
ciency  in  the  cont-ol  design  due  to  43.2%  control  power  wasted  (sq=43.2%  in  Table 
2.2  for  n=8 ,  m=2)  to  truncated  modes  from  an  input  viewpoint.  This  inefficiency 
can  hardly  be  ignored.  One  would  also  want  the  control  design  to  be  efficient 
in  its  control  power,  therefore  assessment  of  spillover  effects  based  on  response 
alone  without  considering  the  control  powers  would  be  premature . 

Fig.  2.4c  shows  the  response  of  the  evaluation  model  both  for  natural 
control  with  continuously  . istributed  inputs  and  LQR  control  with  m-6  point 
inputs.  Both  responses  are  identical.  Because  there  is  no  mode  truncation  in 


29 


the  evaluation  model  the  relative  model  efficiency  is  100%.  On  the  otner  hand, 
the  global  efficiency  is  about  32%  reflecting  the  fact  that  the  LQR  solution  with 
6  point  input  s  uses  about  1  times  more  control  power  than  i t  one  weie  t  o  use  a 
spatially  continuous  input  protile  designed  tor  independent.  eontioL  ot  all  12. 
inodes  In  spite  of  the  more  control  power,  the  LQR  solution  cannot  produce  a 
reap  rse  better  than  that  of  natural  control.  In  the  LQR  solution  with  6  point 
inputs  for  12  modes  intermodal  coupling  of  the  controlled  responses  is  inevi¬ 
table.  In  this  case,  it  is  this  coupling  of  the  controlled  modal  responses  that 
causes  excessive  use  of  control  effort  without  producing  an  improvement  in  the 
controlled  response  over  that  of  natural  control.  This  truly  reflects  the 
inefficiency  of  the  control  design  model. 

Efficient  Model  Reduction 

An  efficient  model  reduction  concept  would  truncate  the  modes  to  which  the 
model  or  global  efficiency  is  least  sensitive.  We  shall  demonstrate  the  model 
redue  ion  technique  based  on  the  model  efficiency.  A  similar  procedure  can  be 
based  on  the  global  efficiency.  However,  for  brevity  we  do  not  demonstrate  this 
alternate  approach. 

We  use  the  model  efficiency  curves  in  Fig.  2.2b  to  find  an  eighth  order 
reduced  model  with  two  and  four  point  inputs.  For  two  inputs  (ni=2)  from  Fig. 
".2b,  we  note  that  the  smallest  increments  in  the  model  efficiency  are  caused  by 
Modes  1,2,  7  and  8.  Hence  we  retain  Modes  3-6,  9-12  as  the  control  design  model 
for  the  given  two  inputs.  Similarly,  for  the  four  input  configuration  (m-4)  from 
tig.  2.2b,  we  note  that  Modes  1,  2,  11  and  12  have  the  least  contributions  to  the 
model  efficiency.  Hence,  we  truncate  these  modes  and  retain  Modes  3-10.  The 
efficiencies  and  control  powers  of  the  new  eighth  order  control  design  models  for 
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two  and  four  inputs  are  shown  in  Table  2.3.  A  comparison  of  these  results  to  the 
control  design  model,  which  was  based  on  the  lowest  eight  structural  modes  (Table 
2.2),  shows  that  the  new  control  design  models  have  significantly  better  effi¬ 
ciencies  and  the  effect  of  model  truncation  for  the  new  reduced  modes  are  insig¬ 
nificant.  For  example,  for  (m=2)  although  the  natural  control  power  S*  has 
increased  from  8  to  24  due  to  the  mode  selection  based  on  model  efficiency,  the 
total  actual  control  power  SR  remained  almost  the  same  (158  vs.  162),  but  the 
modal  control  power  S”  rose  from  90  to  157  which  indicates  that  the  new  control 
design  model  absorbs  almost  all  of  the  actual  control  effort  yielding  a  97%  model 
efficiency.  The  response  profiles  for  the  new  eighth  order  control  design  models 
are  shown  in  Figs.  2.5a  and  b.  Again,  in  these  figures  responses  of  the  corres¬ 
ponding  natural  control,  suboptimal  control,  and  evaluation  model  are  superposed. 
They  are  hardly  different  from  each  other,  the  responses  of  the  suboptimal  con¬ 
trol  and  the  evaluation  model  had  almost  undetectable  overshoots  at  the  peaks  in 
comparison  to  natural  control.  Therefore,  the  curves  were  not  labeled  and  only 
the  response  of  the  evaluation  models  are  shown  in  Figs.  2.5a  and  b.  Among  all 
responses  natural  control  always  achieved  lower  amplitudes  than  the  others. 
Finally,  one  can  now  compare  the  responses  of  the  8th  order  reduced  order  models 
with  m-2 ,  4  to  the  response  of  the  12th  order  evaluation  model  in  Fig.  2.4c  with 
m-6  inputs . 
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Table  2.1:  Efficiencies  for  ACOSS-4  with  n=2 


m 

e*% 

e% 

sq* 

sq% 

1 

0.647792 

1.337612 

2.06488 

152.3056 

98.66239 

2 

1.271 105 

1.271092 

0.99999 

77.67170 

98.72891 

3 

1 . 8iou/l 

1.859323 

1.02269 

53.98065 

98.10468 

4 

2.971592 

2.961611 

0.99664 

32.65535 

97.03839 

5 

3.606529 

3.392659 

0.99615 

26.73133 

96.40734 

6 

4.642842 

4.642761 

0.99998 

20.53855 

95.35720 

Table 

2.2:  Efficiencies  for 

ACOSS-4  with 

n=8 

m 

e*% 

e% 

M 

sq* 

sq% 

1 

2.70735 

56.92516 

21.0262 

15.91033 

43.07480 

2 

5.05710 

56.82256 

11.2362 

8.53799 

43.17744 

3 

19.45213 

63.55967 

3.2675 

1.87333 

36.44030 

4 

29.19059 

69.40089 

2.3775 

1.04825 

30.59911 

5 

39.41595 

69.00026 

1.7506 

0.78648 

30.99740 

6 

47.56904 

72.06676 

1.5150 

0.58721 

27.93324 

Table  2.3:  Effic iencics  and  Control  costs  for  ACOSS-4  with  an  8th  order  model 
obtained  via  Efficient  Model  Reduction  Approach 

Control 

Design 


Mode  s 

m 

e*% 

e% 

4 

sq* 

sq% 

SR 

qH 

cM 

bU 

s* 

3-6,  9-12 

A 

£ 

14 . 7 

96.55 

6.57 

0.24 

3.45 

162.07 

156.47 

5.6 

23.82 

3-10 

4 

44. 7 

93.93 

2.10 

0.14 

6.07 

60.67 

56.99 

3.68 

27.09 
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Figure  2.4b: 


Responses  for  n-8  of  natural  control  (continuous  Input),  and 
suboptlraal  control  and  evaluation  model  with  m-2 ,  e*-5.06%, 
e-56.8%  p-11.2,  sq*-8.3%,  sq-43.2% 
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3.0  EFFICIENCY  MODES  ANALYSIS  OF  STRUCTURE  CONTROL  SYSTEMS 

3.1  Introduction 

In  this  chapter,  we  delve  into  the  structure  of  the  efficiency  quotients 
to  gain  further  insight  to  the  structure -control  problem.  The  developments  of 
this  chapter  center  around  the  fact  that  the  efficiency  quotients  presented  in 
Chapter  2  are  the  ratios  of  quadratic  functionals  of  positive  definite  symmetric 
control  power  functions.  Therefore,  each  one  of  the  quantities  e,  sq,  e*  and  sq* 
represents  a  Rayleigh's  quotient.  It  is  warranted,  then,  to  exploit  the  features 
of  a  Rayleigh's  quotient  to  perform  an  efficiency  modal  analysis  for  the 
structure-control  system  and  study  its  implications. 

In  Section  3.2,  we  discuss  the  computation  of  power  matrices  for  efficiency 
quotients.  Section  3.3  is  central  to  this  chapter.  The  efficiency  eigenvalue 
problem  is  studied  in  this  section,  a  number  of  spectral  properties  of  the  con¬ 
trol  power  matrices  are  discovered,  and  an  efficiency  state-space  transformation 
is  defined.  In  Section  3.4,  the  definitions  of  "controller  modes"  and  "Principal 
Controller  directions"  are  given.  In  addition,  an  efficiency  ellipsoid  is  iden 
tified,  the  surface  of  which  represents  all  initial  disturbances  which  result  in 
a  given  efficiency  of  the  system.  Section  3.5  introduces  the  concept  of  prin¬ 
cipal  efficiency  components  and  discusses  its  role  to  quantify  the  structure  of 
a  link  between  the  controller  design  and  initial  conditions.  Illustrative 
examples,  again  using  the  ACOSS-4  tetrahedral  structure  studies  in  Chapter  2,  are 
given  in  Section  3.6. 

3.2  Computation  of  Control  Power  Matrices  for  Efficiency  Quotients 

Consider  the  N-th  order  evaluation  model  of  the  structural  dynamic  system 


given  by  Eq .  (2.31) 


Mq  +  Kq  -  Q 


(3.1) 


where  q  and  Q  are  the  generalized  coordinates  and  generalized  input  vectors, 
respectively.  Noting  the  transformation  from  a  physical  control  inputs  vector 
F  to  the  generalized  inputs  Q  as 
Q  -  DF 

and  the  modal  properties  given  by  Eqs .  (2.32)  for  the  system  (3.1),  Eq.  (3.1)  can 
be  replaced  by  the  modal  system 

Cn  +  [“2]n£n  -  BnF  .  Bn  -  EjD  (3.2) 

where  £  is  the  N-dimensional  evaluation  model  displacements  vector  of  Eq .  (3.1). 
Typically,  we  shall  assume  that  Eqs.  (3.1)  represents  a  FEM  evaluation  model  of 
the  structural  system.  Introducing  the  2n- dimensional  modal  states  vector  for 
a  control  design  model,  2n  <  2N: 

^  -  [£l  kl  iz  i2---Cn  ^nlT  (3-3) 

we  have  the  modal- state  representation  of  Eqs  (3.1)  and  (3.2)  for  any  smaller 
order  control  design  model 

x  -  Ax  +  BF  (3.4) 

A  —  block-diag  [  ^  ],  B  =  [0  bx  0  b2  0]T  r=l , 2 . . . n<N  (3.5) 

~  cjr  u 

where  br  is  the  r-th  row  vector  of  BN.  The  control  powers  SR,  S”  and  S*  given  by 
Eqs.  ( (2 . 37)  -  (2 .42) ) ,  restated  here  for  ready  reference, 


SR  =  f FTRRFdt  ,  Rr  -  DtM_1D  (3.6) 

S«  «  /FTR«Fdt  ,  R»  =  DTEnE*D  (3.7) 

S*  -  J F*TR*FTdt  ,  R*  =  D*TM_1D*  (3.8) 

become,  for  a  stable  closed-loop  system, 

SR  -  XqPrx0  ,  s£  -  XoPqX0  ,  S*  -  XqP*Xo  (3.9) 

for  which,  the  state  feedback  control  input  is: 
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F  -  Gx 


,  F*  -  G*x 


(3.10) 


The  gain  matrix  G  is  the  control  gain  matrix  available  from  any  control  design 
method.  The  gain  matrix  G*  is  the  gain  matrix  of  the  globally  optimal  solution 
given  by  Eqs .  (2.42),  (2.19)  and  (2.20)  as  discussed  in  Section  2.6. 

The  control  power  matrices  P*  and  P£-  are  obtained  by  solving  the  Lyapunov 
Equations 

Ac[Pr  +  PrAcL  +  GtRrG  -  0  (3.11) 

AjP11  +  P»AcL  +  GtR£g  -  0  (3.12) 

Acl  -  A  +  BG  (3.13) 

and  a  closed- form  solution  for  P*  exists,  as  given  in  Appendix  B. 

Noting  Eq.  (2.38)  and  the  forms  ((3.9),  (3.11)  and  (3.12)),  we  also  have 
for  the  truncated  modes 

S{J  “  ^(P8  -  Pg)  Xq  -  x^PfJxo  (3.14) 

Hence,  the  control  power  wasted  on  the  uncontrolled  dynamics  is  readily  available 
following  solutions  of  Eqs.  (3.11)  and  (3.12)  which  are  based  on  the  controlled 
dynamics  alone,  therefore,  no  knowledge  of  truncated  dynamics  is  required. 

It  must  be  noted  that  the  weighting  matrices  R  in  Eqs.  (3.6)-(3.8)  and, 
hence,  the  control  powers,  are  uniquely  defined  for  computation  of  structural- 
control  system  efficiency.  Having  obtained  the  power  matrices,  one  can  now 
compute  e,  e*,  sq  and  sq*  associated  with  any  control  design  according  to  the 
definitions  (2 . 43) - (2 . 46) . 

3 . 3  Efficiency  Eigenvalue  Problem 

From  the  definitions  of  e ,  sq,  e*  and  sq*,  one  notes  that  they  are  non- 
dimensional  quotients  involving  positive  definite  numerator  and  denominator  power 
matrices  and  the  initial  state  disturbance  Xq.  Therefore,  an  efficiency  quotient 
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is  also  a  Rayleigh's  quotient,  and  an  efficiency  eigenvalue  problem  can  be 
defined  to  analyze  the  structure-control  system. 

In  the  following,  we  only  need  to  define  a  numerator  and  a  denominator 
power  matrix  without  necessarily  referring  to  e,  e*  sq*  and  sq  specifically  and 
use  e  as  a  generic  symbol  to  represent  any  one  of  them.  This  should  not  present 
any  difficulty  as  the  context  will  make  it  clear  whether  we  are  dealing  with  a 
specific  or  a  general  quotient.  Hence,  consider  the  general  form  of  an  effi¬ 
ciency  quotient 

*oPn*o 

e  .  -  P„,  PD  >  0  (3.15) 

XoPdXo 

where  subscripts  N  and  D  denote  numerator  and  denominator  power  matrices.  For 
a  particular  structure  control  system  design  both  PN  and  PD  are  readily  available 
as  given  in  Section  3.2. 

The  efficiency  of  the  given  structure-control  design  system  is  then  depen¬ 
dent  on  the  nature  of  initial  state  disturbance  Xq.  To  see  how  the  initial  state 
interacts  with  the  particular  system,  we  address  the  eigenvalue  problem  associ¬ 
ated  with  the  efficiency  quotient  (3.15). 

Introducing  the  transformation 

x  -  PD‘*z  (3.16) 

an  efficiency  quotient  (3.15)  becomes 

zlEz0 

e - ,  E  -=  Pd'%PnPd-*  (3.17) 

z£z0 

where  E  is  the  nondimensional  symmetric  positive  definite  efficiency  matrix. 

Next,  we  consider  the  standard  eigenvalue  problem  associated  with  the 
matrix  E 

£Vr  -  Arvr  r-1,2 . 2n  (3.18) 
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where  vr  is  a  real  eigenvector  with  positive  real  eigenvalue  Ar,  and  form  the 
modal  matrix  V  of  £ 

V  -  [vj  v2. . ,v2n]  (3.19) 

The  modal  matrix  V  is  normalized  to  satisfy 

VTV  -  I  ,  Vt£V  -  A  (3.20) 

where 

A  -  diag  [Ax  A2...A2n]  (3.21) 

and  I  is  a  2nx2n  identity  matrix. 

Introducing  a  second  transformation 

z  -  Ve  (3.22) 

substitute  it  into  Eq.  (3.17)  and  use  the  normality  conditions  (3.20)  to  obtain 
the  quotient  e  in  the  form 

€o^eo 

e  -  -  (3.23) 

«o«o 

We  shall  refer  to  the  new  states  €  in  Eq.  (3.23)  as  the  efficiency  states.  Com¬ 
bining  Eqs  (3.16)  and  (3.22)  we  have  the  transformation 

x  -  Te  ,  T  -  PD'*V  (3.24) 

In  terms  of  the  efficiency  states,  the  quoti  e  (3.23)  can  be  expanded  in  the 
form 

2n  «oi 

e  -  2  czi\l  ,  cf - ,  0  <  cd  <  1  (3.25) 

i=l  eoeo 

in  which  we  shall  refer  to  c\  and  Xl  as  the  i-th  efficiency  coefficient  and  the 
i-th  characteristic  efficiency,  respectively. 

Next,  introducing  che  inverse  transformation  of  (3.24) 

e  -  Lx  L  -  r1  (3.26) 
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into  Eq .  (3.23)  and  considering  the  original  form  of  efficiency  quotient  in  terms 
of  modal  states  x,  Eq .  (3.15),  we  obtain  the  following  relationships 

PN  -  LtAL  ,  PD  -  LtL  (3.27) 

The  modal-state  dynamics,  Eq.  (3.4)  can  also  be  represented  in  terms  of  the  effi¬ 


ciency  states  by  using  the  transformation  (3.24) 

€  -  Ae«  +  BeF  (3.28) 

A8  =  LAT  ,  B8  =  LB  (3.29) 

The  feedback  control  law  for  F  in  Eq.  (3.10)  becomes 

F  -  Ge  e  ,  Ge  -  GT  (3.30) 

In  terms  of  £-states  the  closed-loop  system  becomes 

€  -  A£L£  ,  A|l  -  A8  +  BeGe  =  LA^T  (3.31) 

where  AqL  is  given  by  Eq.  (3.13). 


By  using  Eqs .  (3.27)  and  (3.31)  in  Eqs .  (3.11)  and  (3.12),  the  Lyapunov 
Equations  associated  with  the  e- state  space  representation  are  obtained  as: 

Acl  A  +  A  A|l  +  G8TR£G8  =  0  (3.32) 

Acl  1  +  A  IgL  +  G8TRrG8  =  0  (3.33) 

Hence,  A  is  recognized  as  the  numerator  power  matrix  and  the  identity  matrix  I 
is  the  denominator  power  matrix  associated  with  the  definition  of  an  efficiency 
quotient  in  the  e -state  space.  That  is 

XoPjjXo  €oTTPn"T€o  foAeo 

e - -  -  -  -  (3.34) 

€oTtPdT€0  £oco 

From  Eq.  (3.34)  or  inverting  Eq.  (3.27)  we  obtain  the  orthonormality  relation¬ 
ships 

TtPnT  -  A  ,  TtPdT  -  I  (3.35) 

Thus,  the  transformation  T  diagonalizes  both  the  numerator  and  denominator  power 
matrices  simultaneously  and  matrix  T  is  orthonormal  with  respect  to  the  real 
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power  control  matrix  P0.  This  completes  the  eigenvalue  problem  of  the  efficiency 
quotient  (3.15) . 


3 . 4  Efficiency  Modes  and  Principal  Controller  Directions  of  a  Structure- 
Control  System 

In  Section  (3.2)  we  showed  that  the  transformation  matrix  T  defined  by  Eq. 
(3.24)  can  be  used  to  transform  from  the  modal  state  space  (x-space)  to  a  new 
efficiency  state-space  (c-space)  to  describe  the  system  dynamics  by  Eq.  (3.28). 
Furthermore,  just  as  the  structural  modal  matrix  E„  diagonalizes  the  structural 
mass  and  stiffness  matrices  according  to  Eq.  (2.32),  the  matrix  T  diagonalizes 
similarly  both  of  the  control  power  matrices  PN  and  PD  according  to  Eq.  (3.35). 
By  analogy,  we  shall  refer  to  the  T  matrix  as  the  efficiency  modal  matrix,  and 
to  its  columns  tr  as  the  efficiency  modes  of  the  control  system. 

T  -  [tx  t2.  .  .t2n]  (3.35) 
Next,  recognizing  the  efficiency  coordinates  cL  in  Eq.  (3.25)  as  the  normalized 
efficiency  coordinates  associated  with  an  initial  disturbance  e0  on  the 
efficiency  state 

Ci  "  «io/|e0|  (3.36) 
and  introducing 

-  Lfx  (3.37) 
where  Li  is  the  i-th  row  vector  of  L-T'1,  the  efficiency  Equation  (3.25)  can  be 
written  in  the  form 


2n  Ai  L^Xq 

E  c? - 1  ,  Ci - 

i-1  e  l*o! 


(3.38) 


Equation  (3.38)  represents  the  equation  of  an  ellipsoid  in  the  2n-dimen- 
sional  space  with  principal  axes  of  length  Ve/Ai  (i-l,...2n)  and  ct  is  the  coor- 
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dinate  along  the  i-th  principal  axis  of  all  initial  disturbances  which  yield  a 
specific  efficiency  e.  Furthermore,  the  direction  of  the  i-th  principal  axis  is 
given  bv  the  i-th  eigenvector  \\  of  the  efficiency  matrix  K  where  V,'s  are 
orthonormal  in  the  usual  sense  as  given  by  (3.20)  associated  with  the  eigenvalue 
problem  of  Eq .  (3.18).  With  respect  to  the  efficiency  eigenvalue  problem  associ¬ 
ated  with  the  original  modal  coordinates  x,  Eq.  (3.15),  the  i-th  principal  axis 
of  the  ellipsoid  (3.38)  is  given  by  the  i-th  column  tA  of  the  transformation 
matrix  T,  Eq.  (3.24)  which  is  the  modal  matrix  of  quotient  (3.15).  We  note  that 
the  directions  tt  are  orthonormal  not  in  the  usual  sense  but  with  respect  to  the 
denominator  control  power  matrix  PD,  given  in  Eq.  (3.35). 

Since  0  <  e  <  1  (excluding  the  case  when  e=sq*) ,  different  initial  con¬ 
ditions  x0  yielding  different  efficiencies  simply  will  result  in  rescaling  of  the 
lengths  of  the  principal  axes  of  the  ellipsoid,  largest  length  in  the  i-th  direc¬ 
tion  given  by  /l/Xi.  The  largest  and  the  smallest  possible  length  of  the  axes 
are  then  /l/Amin  and  /l/Amax,  respectively. 

We  shall  refer  to  the  ellipsoid  (3.38)  as  the  efficiency  ellipsoid,  alter¬ 
nately,  recognize  the  efficiency  modal  vectors  ta  as  the  principal  controller 
directions .  We  note  that  the  efficiency  modes,  characteristic  efficiencies  and 
the  shape  of  the  e-F^f ellipsoid  ?:  completely  and  only  defined  by  the 
particular  control  system  design  and  the  structural  system  design  embedded  in  the 
matrices  G,  A  and  B  in  Eqs .  (3.4)  and  (3.10).  Initial  disturbances  x0  then 
determine  the  size  of  the  efficiency  ellipsoid  acting  merely  as  a  scaling  factor 
to  yield  a  specific  efficiency  for  the  structure-control  system. 

Because  a  quotient  e  in  the  form  of  Eq .  (3.15)  represents  a  Rayleigh’s  quo¬ 
tient  and  since  0  <  e  <  1,  we  observe  the  following  properties  for  efficiencies: 
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a)  The  characteristic  efficiencies  At  are  bounded  by  0  <  Xi  <  1  for  e, 
e*  and  sq,  for  the  quotient  sq*  Ai  >  0. 
b>  Efficiencies  have  stationery  values  at  \t  for  initial  ilistmh.uu.es 
Xq  -  tr  r-1 , 2 , . . . , 2n.  Specifically,  the  minimum  et  t  ic Leucy  a  stLue- 
ture-control  system  can  achieve  is  emin  —  Amin,  an<3  will  occur  if  Xq  - 
tmln,  that  is  if  the  initial  disturbance  coincides  with  the 
eigenvector  (controller  efficiency  mode)  associated  with  Arain,  that  is 
if  the  xQ  is  completely  aligned  with  the  direction  of  the  principal 
axis  along  tmln.  We  shall  refer  to  Amln  as  the  fundamental  efficiency. 
Similarly,  the  upper  bound  of  the  efficiency  of  the  system  is  given 
by  eraax=Amax,  corresponding  to  Xq  -  tmax. 

In  particular,  if  the  quotient  (3.15)  is  the  relative  model  spillover 
quotient 

SyM  XqP^Xq 

sq%  -  100*  -  -  100*  -  (3.39) 

SR  xJp^q 

it  then  represents  the  percent  inefficiency  (fraction  of  control  power  wasted  on 
residual  dynamics)  of  the  control  system.  Noting  that 

PR  -  P£  +  Pjj  (3.40) 

and  denoting  the  eigenvalues  of  the  relative  efficiency  quotient 

XqPcXo 

e  - -  (3.41) 

VT  pRv 
x0r  x0 

by  A®  and  the  eigenvalues  of  the  spillover  quotient  (3.39)  by  Xs  we  write  the 
eigenvalue  problem 

[  A e  PR  -  P«]tr  *  0  r-1,2 . 2n  (3.42) 

substituting  (3.40)  into  (3.42)  we  obtain 
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[  ( 1  -  A®)  PR  -  Pt3]tr  -  0 


(3.43) 


which  constitutes  the  eigenvalue  problem  for  the  relative  model  spillover  quo¬ 
tient  sq .  Hence,  the  eigenvalues  of  the  spillover  quotient  sq  (or  equivalently 
phrased,  controller  power  inefficiency)  are 

A|  —  1  -  AJ  r-1, ,2, . . . ,2n  (3.44) 

with  the  same  eigenvectors  tr  as  that  of  the  relative  model  efficiency  e  in  Eq. 
(3.42).  Similar  to  Eq.  (3.23),  for  the  model  spillover  quotient  we  can  write  the 
expansion 

2n  2n  2n  2n 

sq  -  S  cf  A?  -  Sc?  (1- At)  -  S  c?  -  S  c?  A?  (3.45) 

i-1  i-1  i-1  i-1 

2n 

which  yields,  after  recognizing  S  c?  =  1 

i-1 


sq  -  1  -  e  (3.46) 

as  was  given  in  Eq .  (2.46). 

Note  that  although  there  may  be  infinitely  many  truncated  modes  in  the 
system,  what  happens  from  the  point  of  view  of  control  power  used  is  described 
completely  in  the  2n-dimensional  space  spanned  by  the  controller  efficiency  modes 
tr  r— 1 , 2 . 2n,  for  the  infinite  dimensional  system. 

3 . 5  Principal  Efficiency  Components  -  A  Link  Between  Controller  Design  and 

Initial  Conditions 

The  ultimate  efficiency  of  a  control  system  is  dictated  together  with  the 
structure-control  design  embedded  in  the  matrices  A,  B  and  G  and  the  initial 
disturbance  x0.  To  see  how  these  two  factors  interact,  consider  the  expansion 
of  an  efficiency  quotient  given  by  (3.25)  and  write 
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2n 

e  -  E  ex  .  et  -  cfX,  (3.47) 

i-1 

where  el  represents  the  efficiency  of  the  controller  in  the  i-th  principal  con¬ 
troller  direction.  We  shall  refer  to  e,  as  the  i-th  principal  efficiency  com¬ 
ponent  .  Obviously,  from  Eq .  (3.47),  the  principal  efficiency  ei  is  the  product 
of  the  corresponding  characteristic  efficiency  Aa  and  the  coefficient  cf.  For 
a  given  structure,  the  characteristic  efficiency  Aa  is,  however  only  a  function 
of  the  control  gain  matrix  G  and  the  controller  configuration  represented  in 
matrix  B  via  Eqs.  (3.15),  (3.11),  (3.12)  and  (3.13).  On  the  other  hand,  the 
efficiency  coordinate  cit  Eq .  (3.38)  is  a  function  of  both  the  initial  state  Xq 
and  the  controller  design  via  L  -  T'1  which  in  turn  is  a  function  of  only  G  and 
B,  again  by  virtue  of  Eqs.  (3.15)  and  (3.11),  (3.12)  and  (3.13).  Expressing  the 
functional  dependencies  of  the  quantities  involved  in  the  principal  efficiency 
components  in  (3.47),  one  can  write 

x0Li(A,B,G)Li(A,B,G)x0 

e4  -  cjAA  =  -  Aj.  (A ,  B ,  G)  (3.48) 

Lt(A,B,G)  L(A,B,G) 

and  observe  the  effect  of  both  structure  and  control  design  on  both  Cj  and  Aa. 
Common  sense  requires  that  an  efficient  controller  must  have  high  principal 
efficiency  components  in  the  principal  controller  directions,  in  other  words,  the 
products  cfAj  must  be  high  in  each  direction.  Although  a  control  design  might 
have  a  large  characteristic  efficiencies,  it  will  not  necessarily  have  a  high 
efficiency  unless  the  projection(s)  c4  of  the  initial  disturbance  x0  is  (are) 
significant  along  the  directions  of  large  At .  Considering  that  both  Cj  and  Aj 
are  functions  of  control  design  as  shown  by  Eq .  (3.48),  it  may,  therefore,  be  the 
control  design  itself  which  will  either  enhance  or  degrade  its  efficiency. 
Hence,  it  is  of  significance  not  only  to  have  controllers  with  large 


54 


characteristic  efficiencies  but  also  to  have  their  principal  (Aa)  directions 
aligned  favorably  with  the  initial  disturbances  xG  thus  extracting  large 
projections  of  power  so  that  they  will  yield  high  efficiencies  in  those 
directions . 

Traditionally,  in  the  design  of  a  control  system  (such  as  in  the  LQG-type 
designs)  there  is  no  avenue  to  bring  in  the  initial  disturbance  information  to 
the  computation  of  the  control  gain  matrices.  After  having  obtained  a  so  called 
"optimal"  controller  via  LQG  algorithms,  one  does  not  really  know  how  the  control 
system  will  interact  with  any  initial  disturbance  until  a  simulation  is  done. 
It  is  possible  that  an  optimal  controller  (in  the  sense  of  the  theory  used  to 
design  it)  will  be  an  extremely  inefficient  controller  if  it  does  not  "see"  or 
"receive"  the  initial  disturbance  properly  along  its  principal  controller 
directions . 

On  the  other  hand,  the  efficiency  modes  analysis  of  a  structure-control 
system  clearly  reveals  the  internal  link  by  which  the  controller-design  and  the 
initial  disturbances  xQ  interact.  In  fact,  after  identifying  the  modal  matrix 
T  and  A1,  that  is  the  principal  directions  and  the  characteristic  efficiencies 
of  a  controller  one  can  readily  for  any  initial  disturbance  x0  form  the  projec¬ 
tions  ci  via  Eq .  (3.38)  to  examine  how  they  pair  up  with  the  respective  charac¬ 
teristic  efficiencies  Ai  and  obtain  an  apriori  (before  simulation)  information 
about  the  controller  performance.  Certainly,  the  observation  of  the  link  between 
the  controller  design  and  the  initial  disturbance  is  a  constructive  one  so  that 
given  an  initial  disturbance  and  a  structure,  the  objective  would  be  to  design 
a  control  gain  which  will  pair  up  significant  with  significant  projections  ct . 
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3 . 6  Illustrative  Examples  on  ACOSS-4 


The  structure  is  the  ACOSS-4  structure  considered  in  Chapter  2,  Fig.  2.1. 
In  this  case,  the  structure  is  subjected  to  an  initial  disturbance  of  unit 
displacement  in  the  x  direction  at  Node  2.  Input  configurations  are  selected 
from  a  set  of  twelve  available  actuators  located  at  the  pods  of  the  structure. 
In  all  cases,  the  controls  were  designed  by  using  LQR  theory  with  unit  control 
weighting  (r=l)  and  state  weighting  q=u^  for  the  modal  displacements  and  q-1  for 
the  modal  velocities.  The  evaluation  model  has  N=12  modes  yielding  a  24th  order 
state-space  evaluation  model. 

Efficiency  Modes  Analysis 

Example  3.1:  The  control  design  model  has  two  structural  modes,  Modes  11  and 
12,  n=2 (Modes  11,12)  one  input  was  used  m=l  (#4)  which  was  actua¬ 
tor  number  4.  Figures  labeled  SIM39  (Figures  3.1  through  3.3) 
give  the  line  of  sight  error  simulation  results  for  this  case. 

Example  3.2  n=2  ( Modes :  1 1 , 12 )  ,  m=2(#3,#4),  that  is  two  actuators  #3  and  #4 

were  used  for  control. 

The  nomenclature  for  the  computer  outputs  of  the  examples  presented  in  this 
section  and  the  computer  outputs  are  given  in  Appendix  C.  Simulations  of  the 
globally  optimal  distributed  input  system,  labeled  "Globally  Optimal  System”, 
"Suboptimal  system"  with  m-point  inputs  and  the  "evaluation  model"  response  with 
m-point  inputs  are  presented  for  the  examples  listed  below.  A  designation  SIM 
(number)  associated  with  an  example  indicate  that  figures  labeled  with  SIM  9 
Number  show  the  line  of  sight  error  for  that  example.  Line  of  sight  error  is 
computed  at  the  vertex  (Node  1)  of  the  structure,  according  to  the  expression 
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LOS«U|'1+l,?1 ,  where  Uxl,  Uul  are  nodal  displacements  in  the  x  and  y  directions  at 


the  vertex. 

Example  3.3: 
Example  3.4: 

Example  3 . 5 : 

Example  3.6 


Effect  of  In 
Example  3.7 


n=2  { Modes : 11 , 12 } ,  in“4(#l-4),  SIM41  (Figures  3.4  and  3.6). 

n=8  {Modes:l-8  natural  order},  m-2  (#1,2),  SIM1  (Figures  3.7,  3.8 

and  3.9). 

n=8  (Modes:l-d  natural  order),  m-4  (#1-4),  SIM10  (Figures  3.10, 
3.11  and  3.12) . 

n=2  (Modes: 1-8  natural  order),  m=2  (#1,2).  this  case  presents  the 
results  of  the  (in)  efficiency  modal  analysis  of  the  quotients 
with  PN=Pft,  SIM1  (Figures  3.7,  3.8,  and  3.9). 

itial  Disturbances  on  Efficiencies 

Consider  the  characteristic  efficiencies  for  the  global  efficiency 
e*  and  the  relative  model  efficiency  e  of  Example  3.4  (n»8(l-8), 
m=2)  listed  in  Appendix  C.  For  e*,  we  observe  that  with  a  mere 
change  of  xQ ,  this  system  can  have  a  maximum  global  efficiency  of 
A*ax=53.76%,  corresponding  to  the  2nd  characteristic  efficiency 
listed.  This  would  occur  when  x0-t2,  that  is  when  the  initial 
state  is  the  second  column  of  the  efficiency  modal  matrix  T 
listed.  The  worst  efficiency  will  occur  when  x0-=t1A  corresponding 
to  a  characteristic  efficiency  of  A*in«A^«0 . 773%  as  listed  in  the 
computer  outputs  in  Appendix  C.  For  the  model  efficiency  we  note 
that  all  characteristic  efficiencies  are  59 . 79%>A>54 . 04% .  Hence, 
again  the  corresponding  initial  disturbances  which  will  culminate 
in  respective  efficiencies  can  be  identified  from  the  columns  of 
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the  modal  matrix  T  listed  in  Appendix  C.  It  is  noted  that  regard¬ 
less  of  xOI  the  model  efficiency  can  not  change  more  that  5.7% 
within  that  bracket.  Therefore,  in  this  case  the  model  efficiency 
is  very  robust  to  changes  in  xQ . 
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Figure  3.2:  Suboptlmai  system  LOS  for  n-2  (modes:  11,12),  m-1  (input  #U) 
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Figure  3.3: 


Evaluation  model  LOS  for  n— 2 
example  3.1 


(modes:  11,12}  m-1  (input  #4) 
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4.0  EFFICIENCY  MODEL/CONTROLLER  REDUCTION 


4 . 1  Introduction 

In  the  efficiency  analysis  of  a  2n-th  order  structure-control  system  design 
model,  we  identified  2n-controller  efficiency  modes  which  led  to  the  efficiency 
state  (e -state)  description  of  the  original  design  model  in  the  structural  modal 
space  (x-state).  The  two  modal  descriptions  are  fundamentally  distinct  in  their 
nature.  As  one  is  all  too  familiar  with  structural  modes  of  the  system  charac¬ 
terized  by  the  x-states,  the  controller  efficiency  modes  characterized  by  the  e- 
states  constitute  a  new  concept  where  the  latter  describe  principal  controller 
directions  signifying  avenues  of  efficient  use  of  available  control  power.  From 
a  structure-control  system  perspective,  how  these  two  sets  of  modes  relate  to 
each  other  becomes  a  measure  of  efficiency  of  the  control  task.  A  certain  degree 
of  misalignment  of  these  controller  modal  directions  and  the  structural  modal 
directions  result  in  controller's  not  "seeing"  the  structure  directions  retained 
in  the  model  properly,  thus  not  channeling  the  control  power  to  the  structure  and 
wasting  it  to  the  unrepresented  structure  ;  •'uncated  dynamics.  Conversely, 
certain  structural  directions  may  be  "seen  strongly  or  in  full  view"  by  the 
principal  controller  directions  leading  to  efficient  channeling  of  control  power 
to  the  structure  in  that  direction. 

For  analysis  and  design  purposes  via  the  efficiency  concept,  the  above 
observations  need  to  be  expressed  in  meaningful  mathematical  terms.  To  this  end, 
one  needs  to  correlate  the  structural  modal  states  x  and  the  efficiency  moda* 
states  t .  Quantifying  such  correlations  can  help  the  designer  to  identify  struc¬ 
tural  modal  states  with  weak  or  strong  contributions  to  a  set  of  critical  prin¬ 
cipal  controller  directions.  Once  this  is  accomplished,  one  can  then  proceed  to 
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delete  structural  modal  states  and/or  controller  efficiency  modal  states  which 
are  weakly  correlated. 

The  interrelationship  between  the  structural  modes  and  controller  modes  can 
be  structured  and  quantified  by  decomposition  of  efficiency  quotients  by  defining 
decomposition  matrices  for  the  numerator  and  denominator  control  power  matrices. 
These  decompositions ,  in  turn,  pave  the  way  to  proposing  relevant  model 
controller  reduction  criteria  based  on  the  efficiency  concept. 

In  Section  4.2  we  introduce  decomposition  matrices  for  the  numerator  and 
denominator  control  power  matrices  by  using  the  spectral  properties  of  the 
matrices  given  in  Chapter  3.  Next  component  efficiencies  associated  with  struc¬ 
tural  modal  states  are  defined  for  the  principal  efficiency  components.  Certain 
features  of  these  components  efficiencies  are  discussed  requiring  the  need  to 
define  coupled  and  decoupled  component  efficiencies.  In  Section  4.3,  the  first 
controller  reduction  method  in  the  efficiency  state  space  is  presented.  Section 
4.3  leads  to  the  subject  of  Section  4.4,  efficiency  filtering.  Section  4.5  pre¬ 
sents  the  methodology  for  efficient  modal  controller  reduction  which  relies  on 
the  decomposition  discussed  in  Section  4.2.  Also,  in  Section  4.5  several  impor¬ 
tant  observations  are  made  to  be  able  to  use  the  efficiency  concepts  properly  in 
the  analysis  and  design  studies.  Finally,  in  Section  4.6,  we  give  illustrative 
examples  using  the  ACOSS-4  structure. 

4 . 2  Decomposition  of  Efficiency  Quotients 

Associated  with  any  of  the  quantities  e,  e*,  sq  and  sq*,  we  define  2n- 
dimens ional  numerator  and  denominator  matrices  [ N )  and  [D]  with  elements 

[NJ  -  Njj  -  (XoLiL^AJjj  i  ,  j— 1 . 2 . 2n  (4.1) 

[D!  -  DlJ  -  (4.2) 
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where  Xq  is  the  ini  rial  disturbance  covariance  matrix 

Xc  ~  XqXa  (4.3) 

Next  by  forming  the  total  sums  N  and  D,  one  can  verify  that 

2n  2n 

N  =  S  X  Ny  -  xJPnXo  (4.5) 

i  i 

2n  2n 

D  -  S  S  Dy  -  XqPdXq  (4.6) 

i  j 

e  =  N/D  e  =  ! e , e* , sq , sq* )  (4.7) 

Each  element  Ny  and  Dy  represents  the  power  contribution  of  j  -  th 
structural  modal  state  Xj  in  the  i-th  principal  controller  direction  ^  for  the 
respective  numerator  and  denominator  terms.  We  shall  refer  to  matrices  [N]  and 
[D]  as  the  numerator  and  denominator  control  power  decomposition  matrices, 
respectively.  The  quotient  e  in  terms  of  sums  of  their  elements  represents  an 
efficiency  decomposition.  Although  through  the  controller  design  one  has  the 
power  matrices  PN  and  PD  available  in  the  beginning,  their  representations  in 
terms  of  the  elements  Ny  and  Dy  unveils  the  internal  structure  of  the  power 
associated  with  the  controller.  This  exposition  of  the  internal  structure  is, 
however,  possible  only  through  the  efficiency  modal  analysis.  Specifically 
through  this  modal  analysis,  each  element  Ny  and  Dy  quantifies  the  correlation 
between  the  j  -  th  structural  mode  Xj  and  the  i-th  controller  mode. 

Next,  we  define  the  partial  sums 

?  n  2  n 

s  Nu  •  Dj  -  X  Dy  j-1,2 . 2n  (4.8) 

i-1  i-1 
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as  the  power  contributions  (extraction)  of  j-th  structural  modal  state  Xj  to 
(from)  all  principal  controller  directions.  Similarly,  we  define  power 
contributions  (extraction)  of  all  structural  modal  states  to  (from)  the  i-th 
principal  controller  direction: 

2n  2n 

Ni-2  Ntj  ,  Dt  -  Z  Du  i-1,2 . 2n  (4.9) 

J-l 

Considering  definitions  (4.9)  and  the  definition  of  principal  efficiency  com¬ 
ponent  eit  Eq.  (3.47),  it  follows  that 

Ni  N, 

ei  -  -  - - c2iXi  i-1,2 . 2n  (4.10) 

ZDi  D 

i 

Where  we  recall  that  eL  represents  a  component  with  respect  to  the  i-th  principal 
controller  direction.  Noting  Eq.  (3.25),  one  can  also  regard  et  as  the  contribu¬ 
tion  to  efficiency  of  the  i-th  efficiency  state  tj.  Similarly,  one  can  identify 
efficiency  components  with  respect  to  the  structural  modal  directions  and  define 

NJ  Nj 

ej  -  -  -  -  ,  e  -  Ze,  j-l,2,...,2n  (4.11) 

D  ZD,  j 

j 

We  shall  refer  to  e,  as  the  j-th  (structural)  component  efficiency  and  emphasize 
the  difference  in  the  terminology  that  is  being  used  here  that  ej  is  the  com¬ 
ponent  efficiency  and  ei  is  the  principal  efficiency  component,  respectively,  the 
former  is  with  respect  to  structural  modal  components  x,  and  the  latter  with 
respect  to  controller  efficiency  modal  components  eL. 

Strictly  speaking,  although  definitions  of  components  ei  and  e,  as  in  Eqs . 
(4.10),  (4.11)  are  mathematically  correct  since  their  sums  over  respective 
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indices  yield  the  total  efficiency  e,  from  a  physical  perspective  certain  fea¬ 
tures  of  these  definitions  must  be  considered  cautiously.  In  both  definitions, 
the  denominator  is  common  for  all  components  e.  and  <=r  .  So  in  effect,  the  sig¬ 
nificance  of  a  component  is  solely  determined  by  the  numerator  power 
contributions  NA  and  Nj .  For  that  matter,  the  use  of  D  which  is  the  real  control 
power  SR  as  a  common  denominator  is  of  no  consequence,  as  the  relative 
contributions  of  Nj  (or  Nt)  terms  to  total  efficiency  would  not  change  if  we  were 
to  use  any  other  quantity  as  a  common  denominator. 

Indeed,  in  this  case  the  function  of  efficiency  as  a  relevant  nondimen- 
sional  quantity  is  subordinated  to  the  contributions  to  the  dimensional  numerator 
control  powers.  Therefore,  from  the  efficiency  point  e^  and  e^  component  des¬ 
criptions  may  not  be  proper  characterizations.  In  case  of  component  efficiencies 
ej  for  the  structural  modal  states  Xj,  the  characterization  is  completely  out  of 
contact  with  the  concept  and  purpose  of  efficiency  quotient  as  we  show  next. 

Considering  the  definition  of  Nj ,  Eqs .  (4.8)  and  (4.1)  and  the  spectral 

decomposition  of  control  powers  matrices  in  (3.27),  we  can  write 

2n 

Nj  -  [Xq  2  LjAjLl] jj  =  [XqPnIjj  =  ejD  j-1,2 - 2n  (4.12) 

i-1 

summing  over  all  component  efficiencies: 

2n 

eD  —  2  Nj  =  Trace  XgPjj  =  XqP^Xq  =  (4.13) 

J 

Thus,  component  efficiency  characterization  ej  in  this  sense  is  no  different  that 
component  cost  analysis  of  a  modal  power  s£  of  the  control  design  model  as 
defined  in  Ret.  (7)  and  all  the  benefit  and  the  physical  reason  in  introducing 
a  denominator  control  power  matrix  to  yield  the  nondimens ional  efficiency  is 
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lost.  Furthermore,  the  components  Nj  given  by  Eq.  (4.12)  are  computed  exactly 
as  given  in  Ref.  (7)  associated  with  which  there  is  certainly  no  concept  of 
efficiency. 

A  similar  situation  exists  for  the  efficiency  components  eit  however,  at 
a  conceptually  less  critical  level.  Considering  the  definition  of  Nif  Eqs . 
(4.9),  (4.2)  and  (3.27),  we  can  write 

2n 

Ni  =  ^(XoLiAiLDjj  i=l , 2 . 2n  (4.14) 

j-1 

which  upon  recognizing  XqLi  -  tlQ  and  E(xoL*)jj  -  «io  yields 

j 

Ni  “  Ajffo  “  eiD . 

Summing  over  i 

2n 

SNi  -  SAief0  -  eD  -  Sg  -  «5a«0  (4.15) 

i  i 

This  is  again  a  decomposition  of  the  numerator  modal  power- of  the  control  design 
model  regardless  of  the  function  of  denominator  control  power  D  as  the  nondimen- 
sionalization  factor.  The  form  of  S”  as  in  Eq .  (4.15)  is  similar  to  the  cost 

decoupling  concept  discussed  in  Refs.  (7  and  12).  Here,  the  efficiency  states 
t  play  a  role  similar  to  the  cost  decoupling  coordinates.  Hence,  principal 
efficiency  components  e^  characterization  of  e  is  also  a  characterization  of  the 
numerator  control  power  S”  alone.  However,  because  this  characterization  can  be 
af tec  led  caly  in  terms  of  efficiency  coordinates  ei  after  solving  for  the 
efficiency  modal  matrix  T,  the  function  of  denominator  control  power  in  the  con¬ 
cept  of  efficiency  is  accounted  for  indirectly.  In  contrast,  in  the  ej  charac¬ 
terization  there  was  no  need  for  the  efficiency  concept.  It  must  be  noted,  how¬ 
ever,  that  the  identification  of  the  cost  decoupling  coordinates  as  in  Ref.  (7 
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and  12)  is  only  a  recognition  of  the  algebraic  similarity.  The  mathematical  and 
the  physical  foundations  under  which  the  decoupling  features  are  brought  about 
in  this  investigation  and  in  Ref.  (7  and  12)  are  fundamentally  different. 

From  the  above  discussion,  to  be  consistent  with  the  definitions  of  effi¬ 
ciency  there  evolves  the  need  to  define  component  efficiencies  e^  for  the  contri¬ 
butions  of  structural  modal  states  Xj  which  takes  into  account  the  denominator 
contributions  Dj  as  well  as  the  numerator  contributions  Nj .  To  this  end,  we 
define  two  component  efficiencies.  The  first  one  is: 

NjD  +  DjN  2n 

6j  - - e  -  2  ej  (4.16) 

2D2  j 

which  we  shall  refer  to  as  the  coupled  component  efficiency.  The  coupled  desig¬ 
nation  is  used  to  note  that  in  the  definitions  of  Nj  and  Dj ,  terms  involving 
products  of  xJO  with  all  other  states  appear.  Hence,  a  state  xjo  may  have  a  high 
component  efficiency  contribution  not  because  of  itself  alone  but  also  because 
of  the  contribution  of  its  coupling  terms  with  another  state  or  states.  The 
second  definition  of  component  efficiency  disregards  all  coupling  terms  among  the 
states,  and  we  refer  to  it  as  the  decoupled  component  efficiency: 

ei  "  nj/dJ 

nj  ”  Xjo^Njjxjo 

(4.17) 

dj  “  Xjo^Dj  jXjo 

XI  -  Uj  £j]  j-1 . n 

where  PNjJ  and  PDjj  are  the  2x2  block  diagonal  partitions  of  the  numerator  and 
denominator  power  matrices,  respectively.  In  the  definition  of  component  effi¬ 
ciency,  we  retained  the  pair  of  a  modal  displacement  and  its  velocity  as  a  com¬ 
ponent  for  obvious  reasons.  This  is  also  to  be  observed  in  all  component  effi- 
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ciencies.  The  decoupled  component  efficiency  ej  would  be  exactly  the  efficiency 
of  the  system  if  all  initial  states  but  the  modal  state;;  j  (displacement  and 
velocity)  were  zero,  hence,  they  indeed  represent  proper  component  contributions 
to  the  overall  efficiency.  We  also  recognize  that  ej  <  1 . 

Next,  noting  that 


N  —  Snj  +  S  SxJ0PNijX0j  i.j-1,2 . n  (4.18) 

j  i*j 

and 

D  “  Sdj  +  S  SxioPDijX0j 

j 


where  the  second  summations  are  the  off  block-diagonal  power  cross-coupling 
terms.  In  general,  these  off-diagonal  contributions  are  sign  variant,  and 
therefore,  one  might  expect  them  to  vanish  in  a  statistical  sense  if  sufficient 
modal  states  are  disturbed  initially.  Thus,  we  can  write  an  expected  total 
efficiency  in  terms  of  the  block-diagonal  decoupled  powers  as 


n 

Zrij 

j 

E(e)  =  -  (4.19) 

n 

Sdj 

j 


where  E  denotes  expectation.  Further,  we  caution  that 


n 

e  *  Eej. 

j-1 


Many  studies  that  have  been  conducted  during  the  course  of  this  investigation 
verify  that  the  expectation  of  overall  efficiency  in  terms  of  decoupled  component 
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contributions  is  quite  relevant  and  off-block  diagonal  contributions  to  both  the 
numerator  and  control  powers  are  insignificant. 

In  the  definitions  of  both  the  coupled  and  decoupled  component  efficiencies 
(ej)  one  notes  that  the  numerator  and  denominator  component  contributions  Nj ,  Dj , 
nj  an  dj  can  be  computed  exactly  directly  from  the  numerator  and  denominator 
power  matrices  PN  and  PD  without  invoking  any  of  the  results  of  the  efficiency 
modes  analysis.  In  these  definitions  of  the  exact  component  contributions  via 
the  power  matrices  summation  over  all  principal  controller  directions  i— 1 , 2 , . . , 2n 
is  implicit.  However,  in  reality  not  all  principal  controller  direction  contri¬ 
butions  are  significant,  that  is  the  least  efficient  controller  directions  need 
not  be  taken  into  account,  therefore,  it  suffices  in  general  to  consider  only  the 
partial  sums  involving  dominant  controller  modes  in  Eqs .  (4.1),  (4.2)  and  (3.27) 

np  np 

Nj  -  [XqS  PHJ]jj  ,  Dj  -  [X0  S  PDi]jj,  np  <  2n  j-1 . 2n(4.20) 

i-i  i-i 

where 

PN1  =  ,  PDi  =  LtLl  (4.21) 

are  controller  power  spectral  components  of  the  power  matrices  available  only 
after  an  efficiency  modal  analysis.  Similarly,  for  approximations  to  ^  and  dj : 

n?  nP 

^Njj  ~  ^  [  Phi  J  j  j  an<*  ^Djj  “  S  [PDi]jj,  j=l,2,...,2n  (4.22) 

i  i 

We  are  now  in  a  position  to  propose  model/controller  reduction  criteria  based  on 
the  concept  of  efficiency.  To  this  end,  we  shall  utilize  the  component  effi¬ 
ciencies  ej  and  principal  efficiency  components  ej  discussed  in  this  section. 


79 


4 . 3  Controller  Reduction  in  the  Efficiency  State  Space  -  Efficiency  Components 

Truncation 

We  use  the  transformation  (3.26)  from  the  modal  space  to  the  efficiency 
space  to  represent  the  closed- loop  system  in  the  efficiency  space  as  given  by 
Eqs .  (3 . 28) - (3 . 31)  . 

A  controller  reduction  in  the  e- space  can  be  affected  by  retaining  only  nP 
components  of  e  associated  with  the  dominant  controller  modes  with  highest  rip 
principal  efficiency  components  eL ,  Eqs .  (3.47)  that  is  { e± } x  —  {e1>e2>  ...  enp)  , 
rip  <  2n.  Rearranging  the  controller  modal  matrix  T  in  accordance  with  the  set 
(ei)1  we  write: 

x  —  T1e1  +  T2«2  (4.23) 
where  Tj  are  the  controller  modes  corresponding  to  the  characteristic  effi¬ 
ciencies  Xi  associated  with  the  set  (e1}1  and  subscript  2  denotes  the  truncated 
states.  The  decreasing  order  of  principal  efficiency  components  (ei)  does  not 
necessarily  correspond  to  highest  rip  characteristic  efficiencies  Xi .  Similarly, 
the  partitioned  form  of  the  closed-loop  system  in  the  efficiency  space  is: 


el 

Afi+BfC?  • 

a|2+b®g2 

«l 

Afi+BfG®  : 

Alz+BlGf 

.*2  . 

(4.24) 


where  the  matrices  A®,  B®  and  G®  are  given  by  Eqs.  (3 . 28)  -  (3 . 31)  .  A  reduced 
order  controlled  system  is  obtained  from  the  above  by  considering  only  the 
dynamics  associated  with 


=  f  A 


li 


^l^i  1  €  i  =  A:lii€  i 


;  10 


=  e 


10 


(4.25) 


ior  which  the  feedback  control  law  is  obtained  by  rendering  Gf  -  0  in  the  full 
order  control  law  Eq .  (3.30). 


F  -  Gfi ! 


(4.26) 
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However,  considering  the  Cull-order  2n- th  order  design  model  Eq .  (4.24)  with 

G2=0 ,  the  actual  system  dynamics  with  the  reduced-order  controller  would  be, 


1 1 

Afi+BJGf 

1 

i  > 

K> 

1 _ i 

£1 

.‘2 . 

Af1+BfGf 

a22 

iz 

(4.27) 


Thus,  in  the  actual  system,  the  reduced-order  states  e!  and  the  truncated  states 
r 2  retain  their  dynamic  coupling.  Because  the  simulations  of  (4.25)  and  (4.27) 
will  yield  different  responses  for  Cj,  an  overbar  is  used  for  ex  in  (4.25)  in 
ignoring  the  dynamic  coupling  of  with  the  truncated  efficiency  states  e2.  The 
implicit  assumption  in  this  approach  is  that  the  system  (4.27)  remain  stable  with 
the  reduced- order  controller.  This  assumption  seems  reasonable  since  the  reduc¬ 
tion  criterion  is  based  on  retaining  the  states  associated  with  the  most  effi¬ 
cient  controller  directions,  thus,  the  anticipation  is  that  significant  control 
power  still  remains  with  the  reduced-order  controller.  Put  in  other  words,  the 
truncated  controller  directions  associated  with  the  e2  states  with  low  principal 
efficiency  components  et  will  have  high  inefficiencies,  thus,  their  elimination 
will  save  control  authority.  On  the  other  hand,  the  stability  of  the  reduced 
model  Eq .  (4.25)  is  guaranteed  as  we  will  note  shortly.  In  addition,  the  effi¬ 
ciency  of  the  reduced  order  model  (4.25)  with  the  reduced-order  controller  (4.26) 
would  be  of  interest.  Assuming  that  the  Aclii  is  stable,  the  Lyapunov  Equations 
associated  with  the  modal  power  of  the  design  model  S”  and  the  real  power  SR 
become 

Aclii  +  P^lh  +  CfR»Gf  «  0  (4.28) 

Aclii  +  P^LU  +  GfTRRGf  -  0  (4.29) 
Next,  considering  the  nPxn  upper  left  hand  partitions  of  Equations  (3.32)  and 
(3.33),  we  obtain 
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Aclii  ^11  +  AnA£L11  +  G®TRcCf  -  0 

Aclu  In  +  IuA^lh  +  GfRRGf  -  0 

Comparing  Eqs .  (4 . 28) - (4 . 31)  we  conclude: 

pM  _  A  pR  _  T 

re  _  A11  •  re  _  ill 

where 


(4.30) 

(4.31) 


(4.32) 


An  =  diag[AJ  A®... A®]  (4.33) 
and  Iu  is  the  np-dimensional  identity  matrix.  Conversely,  associated  with  the 
closed-loop  matrix  A^m,  since  Eqs.  (4.30),  (4.31)  guarantee  An  and  In  as 
positive-definite  solutions,  it  follows  that  AcL11  is  a  stable  matrix.  Thus,  we 
have  shown  that  the  reduced-order  modal  matrix  Tp  still  simultaneously  diagonal¬ 
izes  the  modal  and  real  power  matrices  s£-  and  SR  of  the  reduced- order  model 
(4.25)  without  perturbing  the  original  characteristic  efficiencies  associated 
with  the  states.  The  nP-th  order  efficiency  of  the  reduced  system  becomes 


enP 


e  iqAu€  10 


(  I0e  10 


nP 

-  2  A® 


1-1,2, 


(4.33) 


The  change  in  the  efficiency  of  (4.25)  in  comparison  to  the  original  system 
(4.24)  is 

2n  2n 

Ae  -  e-enP  -  E  c^A®  -  E  et  (4.34) 

i=nP+l  i=nP+l 

Hence,  the  efficiency  of  the  reduced- system  order  in  Eq .  (4.25)  will  decrease  by 
the  amount  of  truncated  efficiency  components  of  the  original  system  of  Eq . 
(4.24)  . 

This  method  of  controller  reduction  in  the  e -space,  because  of  the  diagonal 
forms  of  the  power  matrices  Au  and  I1P  is  akin  to  the  controller  reduction  via 
cost  decoupled  coordinates  discussed  in  Ref.  (7). 
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The  controller  reduction  in  the  efficiency  space  is  attractive  from  the 
point  of  view  of  that  the  reduced-order  controller  (4.26)  requires  no  new  com¬ 
putations  for  the  control  design  since  it  is  readily  obtained  by  disregarding  the 
control  gains  associated  with  the  £2-states  in  original  control  law  (3.30). 
Furthermore,  the  efficiency  of  the  reduced  system  (4.25)  requires  no  new 
computation  either,  since  it  is  readily  given  by  truncating  the  efficiency 
components  of  the  original  system  at  i-np  in  (4.33).  The  method  also  has  some 
disadvantages.  The  reduced  order  system  (4.25)  neglects  dynamic  coupling  with 
the  £2- states  which  has  to  be  assessed  although  one  hopes  that  truncation  of  the 
least  efficient  states  t2  results  in  truncation  of  least  important  dynamic 
coupling  terms.  Because  of  ignoring  this  dynamic  coupling  which  is  referred  to 
as  dynamic  spillover  between  £P  and  e2,  the  sole  assessment  of  the  reduced-order 
controller  (4.26)  can  only  be  obtained  by  simulation  of  the  2n-th  order  system 
(4.27)  instead  of  the  nP-th  order  system.  However,  in  this  case  one  notes  that 
the  reduced  control  to  be  simulated  in  (4.27)  has  the  form  F-G®£x  instead  of 
F-G®^.  If  one  were  to  evaluate  (4.27)  by  using  the  latter  form  of  reduced 
control  using  feedback  of  this  would  constitute  an  open-loop  control  with 
respect  to  the  £L  and  £z  dynamics  and  the  system  (4.27)  would  still  remain  neu¬ 
trally  stable.  The  problem  is  that  although  the  £ -states  (or  any  other 
cost/power  decoupling  states  as  presented  in  Refs.  11,  12)  diagonalize  the  power 
matrices  PN  and  PD  they  remain  dynamically  coupled.  The  truncation  of 
dynamically  coupled  states  always  leads  to  uncertainty  of  the  behavior  of  the 
retained  states  in  the  actual  system.  In  contrast,  truncation  of  dynamically 
decoupled  states  still  insures  that  the  retained  states  in  the  actual  system 
behave  exactly  as  they  do  in  the  reduced-order  model  of  the  system.  Certainly 
structural  modal  coordinates  have  precisely  this  feature  and  this  should  point 
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out  why  l  hey  «-i  it-  ever  so  dcs  i  r  at)  1 o  i  ti  tin*  ana  1  y  s  l  s  aiul  synthesis  of  structural 
dynamic  systems.  This  being  the  case,  it  makes  practical  sense  only  to  consider 
the  reduced  order  controller  form  as  F=G“e1  instead  of  using  feedback  of  tl  in 
(4.25)  which  leads  one  to  consider  Eq .  (4.27)  which  in  turn  constitutes  a  closed- 
loop  system.  Finally,  because  the  efficiency  ep,  Eq.  (4.33)  of  the  reduced-order 
system  given  by  (4.25)  is  obtained  without  regard  to  dynamic  spillover  between 
and  £ 2  states,  the  true  efficiency  of  the  dynamics  of  e1  states  in  (4.27)  when 
the  dynamic  spillover  is  considered  will  be  different  from  eP,  ta  (4.33).  This 
effect  can  indeed  be  significantly  different  in  some  cases  so  that  the  simplistic 
optimism  provided  by  Eq .  <4.34)  to  anticipate  the  change  in  the  efficiency 

becomes  unrealistic. 

The  evaluation  of  the  stability,  efficiency  and  the  response  of  the  system 
(4.27)  will  necessitate  a  recycle  of  the  stability  and  efficiency  analysis.  This 
being  the  case,  it  proves  more  convenient  to  readdress  the  £ -state  truncation 
approach  in  terms  of  the  original  modal  states  x  of  which  the  structure-control 
designer  has  a  better  understanding.  In  this  form,  when  studied  in  the  modal 
space,  the  re-evaluation  of  the  e -space  controller  reduction  is  tantamount  to  an 
efficiency  filtering  method.  This  we  show  in  the  next  section. 

4 . 4  Efficiency  Filtering 

Implementation  of  the  reduced-controller  in  the  form  of  Eq .  (4.26)  requires 
estimation  of  the  e,  states  which  will  theoretically  be  different  than  those 
states  £t  in  (4.27)  with  dynamic  spillover  considered.  However,  from  an  imple¬ 
mentation  point  of  view,  there  is  no  advantage  gained  in  feeding  back  the  £j  as 
discussed  in  the  previous  section,  as  opposed  to  feeding  back  the  actual  el 
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states  with  which  there  is  an  advantage  over  for  the  reasons  just  d'scussed. 
Thus,  ve  consider  the  reduced  feedback  controller 

F  “  G!*i  (4.35> 

Then,  Eq .  (4.27)  represents  the  actual  system  in  the  e-space.  However,  the  con¬ 
trol  law  (4.35)  ana  the  system  (4.27)  can  be  transformed  back  to  the  original 
modal  x-space.  write  ^he  control  law  (4.35)  as 

F  -  [Gf  0]  I  e jj  (4.36) 

and  note 

(Gf  Gf]  -  r[ T,  :2]  (4.37) 

Introducing  the  portions  into  (4.36),  we  obtain 

F  -  GT1L1x  -  Gx  (4.38) 

which  yields  the  new  closed- loop  dynamics 

x  »  [A+BGjx  =  [A+BGTiLJx  (4.39) 

Therefore,  the  e -space  truncation  of  e2  states  is  equivalent  to  full  x-state 
feedback  through  the  new  feedback  gain  matrix  G  which  simply  represents  a  pro¬ 
jection  of  the  original  G  matrix  through  the  matrix  T^  which  i:;  readily  avail¬ 
able  as  soon  as  the  partition  is  decided.  The  matrix  T1L1  essentially  filters 
out  the  contribution  of  the  least  efficient  e2  states  to  the  x-states.  Hence, 
we  refer  to  this  method  as  efficiency  filtering. 

Alternately,  noting  that 

TiLi  "  i'T2L2  (4.40) 

the  closed- loop  system  becomes: 

*  -  [A  *  BG 1 1  -  T2L2)  1  x  -  (Acl  -  BGT2L2  ]  x  (4.41) 

or 

x  -  [Acl  +  AG ] x  ,  AG  -  BGT2L2  (4.42) 
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hence,  the  closed- loop  systems  oi  Eqs .  (4.39)  and  (4.42)  are  equivalent  and  they 
can  also  be  viewed  as  the  perturbation  of  the  original  system  (3.31)  through  the 
modal  gain  error  matrix  AG.  Because  all  matrices  involved  in  AG  are  known,  the 
gain  perturbation  is  a  structured  perturbation.  The  stability  and  the  perform¬ 
ance  robustness  (perturbation  in  P£?  and  PR)  of  the  perturbed  system  v4.42)  can 
now  be  studied  by  the  methods  of  Refs.  (12)  and  (13).  From  the  performance 
bounds  on  Pj3  and  PR,  the  robustness  of  efficiency  to  filtering  the  e2  states  can 
easily  be  inferred.  We  remind  that,  the  study  of  stability  and  performance  of 
efficiency  filtered  closed-loop  system  (4.42)  in  terms  of  the  x-state  is  equiva¬ 
lent  to  the  study  of  the  e -state  Equation  (4.27)  by  using  the  reduced-controller 
(4.35). 

The  « -state  truncation  and  the  efficiency  filtering  technique  as  a  way  cf 
obtaining  and  evaluating  reduced  order  controllers  do  not  seem  attractive  since 
the  approach  anticipates  a  lower  efficiency  in  the  reduced  system.  Furthermore, 
the  truncation  of  efficiency  states  represents  a  further  degree  of  abstraction 
in  physical  understanding  of  the  controller-reduction  in  contrast  to  truncation 
of  modal  states.  Also,  there  is  the  problem  of  dynamic  spillover  associated  with 
the  approach,  therefore,  it  is  plausible  to  search  for  controller  reduction  cri¬ 
teria  which  will  culminate  in  an  increase  in  the  efficiency  of  the  reduced-order 
system  and  the  reduction  is  affected  in  terms  of  the  original  modal  (x-space) 
states.  At  the  very  least,  if  no  significant  degradation  of  the  efficiency  is 
experienced  in  such  a  reduction  process,  the  method  would  be  preferable  over  the 
€ -state  reduction  method.  The  underlying  philosophy  of  such  a  method  is  then  to 
identify  and  truncate  modal  coordinates  x  in  the  design  model  which  have 
significant  inefficiency  associated  with  them  so  that  their  deletion  will  make 
more  efficient  use  of  control  power  available  to  the  modes  retained  thereby 
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resulting  in  an  increase  in  the  efficiency  of  the  controller  despite  a  reduced- 
order  system. 

4 . 5  Efficient  Modal  Controller  Reduction 

The  premise  of  the  controller  reduction  methodology  we  shall  propose  in 
this  section  is  that  the  structural  modal  states  Xj  that  contribute  least  to  a 
given  set  of  principal  efficiency  components  eA  along  the  associated  principal 
controller  directions  are  to  be  deleted.  Conversely,  the  modal  coordinates  which 
contribute  most  to  a  given  set  of  spillover  components  sqi  along  the  associated 
principal  controller  directions  are  to  be  deleted.  It  turns  out  the  two  views 
do  not  necessarily  yield  the  same  set  of  truncated  modal  states. 

The  method  involves  computation  of  coupled  or  decoupled  component  effi¬ 
ciencies  ej  as  defined  by  Eqs .  (4.16)  and  (4.17)  for  an  initial,  structural- 
control  design  model  and  performing  the  2n-th  order  efficiency  modal  analyses  for 
one  or  more  of  the  quotients  e,  e*  sq  and  sq*.  Upon  studying  the  characteristic 
eigenvalues  and  the  principal  efficiency  components  of  these  quotients  along  the 
principal  controller  directions  the  significant  components  of  relevant  efficiency 
and/or  spillover  quotients  are  recognized  to  retain  rip  <  2n  principal  directions 
in  the  computation  of  component  efficiencies  ej  according  to  Eqs.  (4 . 20) - (4 . 22) 
where  the  summation  over  i  involves  rip  controller  directions.  Note  that  depend¬ 
ing  on  which  quotients  are  studied,  the  quotient  ej  in  Eqs.  (4 . 20) - (4 . 22)  could 
specifically  represent  any  of  quantities  e j ,  ej,  sqj  and  sq^f. 

After  computing  the  component  efficiencies,  one  then  retains  a  desired 
number  nR  of  the  structural  components  with  the  highest  component  efficiencies 
ej  and/or  e* .  Conversely,  if  the  above  analyses  are  based  on  spillover  quotients 
one  retains  nR  structural  modal  coordinates  with  the  least  component  spillover 
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quotients  sq£  and  sq*.  As  shown  in  Section  (3.3)  e  and  sq  have  the  same 
eigenvectors,  that  is  the  same  principal  controller  directions  and  their 
eigenvalues  are  related  as  (  \“ )  and  {As}  ■-  (l-X*).  Therefore,  if  e  and  sq  are 
considered,  modal  analysis  ot  either  one  wilL  he  sufficient.  It,  however,  e*  and 
sq*  are  desired,  no  immediate  relationship  between  their  eigenvectors  and 
eigenvalues  are  recognized,  so  modal  analyses  for  both  of  them  must  be  performed. 
Furthermore,  recall  that  since  e,  e*  and  sq  cannot  be  greater  than  unity  their 
characteristic  eigenvalues  cannot  be  greater  than  unity.  The  eigenvalues  of  the 
global  spillover  quotient  sq*  can  be  greater  than  unity.  After  selecting  the  nR 
structural  modal  coordinates  as  suggested,  two  avenues  are  available.  Restating 
here  for  convenience  the  original  control  design  model 

x  -  Ax  +  BF  ,  F  -  Gx 

on  which  tne  efficiency  analyses  are  performed,  the  first  avenue  is  to  retain 
only  the  2nR  states  xR  and  truncate  the  remaining  modal  states  and  implement  the 
reduced  order  model 

xr  =  Arxr  +  BrF  ,  F  =  Grxr  2nR  <2n  (4.43) 
without  redesigning  the  feedback  control  law  where  AR,  BR  and  GR  are  the  parti¬ 
tions  of  A,  B  and  G  corresponding  to  the  2nR  states  xR.  This  is,  in  general,  the 
most  common  spirit  of  controller  reduction.  To  assess  the  effectiveness  of  the 
reduced  order  model  the  stability  and  new  efficiency  of  the  system  (4.43)  must 
be  studied.  In  general,  the  characteristic  efficiencies  and  the  modes  of  the 
controller  will  change.  But,  if  the  reduced-order  model  identified  is  a  good  one 
its  efficiency  features  will  have  changed  over  the  original  system  in  the 
direction  that  the  analysis  anticipates,  that  is  truncation  of  least  efficient 
or  most  inefficient  modal  states  should  leave  behind  at  least  a  more  efficient 
reduced  model. 
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The  second  avenue  available  for  controller  reduction,  rather  than  trun¬ 
cating  the  original  control  law,  is  to  redesign  the  feedback  control  law  for  the 
set  of  modal  states  xR 

*r  =  ar  +  BrFr  ,  Fr  -  GrjjXr  (A./44) 
where  Grr  is  the  new  redesigned  control  gain  matrix  for  the  reduced-order  model 
for  the  reduced  structural  system.  For  example,  if  the  initial  control  design 
is  obtained  via  LQR  it  will  be  repeated  for  the  system  with  the  performance 
measure 

2J  -  /(xJW,*  +  FRWCRFR)dt  (4.45) 
Of  course,  the  control  redesign  is  a  conservative  approach.  It  guarantees  the 
stability  of  the  reduced  order  model  and  its  efficiency  can  be  reanalyzed. 
Intuitively  one  would  stand  a  better  chance  of  improving  it  with  the  control 
redesign. 

The  redesigning  of  the  controls  is  viewed  here  as  a  more  realistic  and,  in 
fact,  practical  approach  to  the  model/controller  reduction.  The  reduced  order 
model  (4.44)  is  smaller  than  the  original  control  design  model  for  which  algor¬ 
ithmic  and  theoretical  capability  exists  to  begin  with.  Consequently,  applica¬ 
tion  of  the  same  capability  to  a  reduced  order  system  by  using  the  same  tools 
should  not  be  a  problem.  Furthermore,  the  process  of  identifying  a  good  reduced 
order  model  is  an  off-line  procedure  and  performing  control  redesigns  and  effi¬ 
ciency  modal  reanalyses  do  not  require  more  computational  resources  than 
available  to  perform  the  original  2n-th  order  design.  Countless,  aimless 
redesigns  are  not  advocated  here  to  identify  a  good  reduced  model  that  satisfies 
the  system  constraints.  On  the  contrary,  the  efficiency  modal  analysis  presented 
heretofore,  does  lay  a  foundation  for  the  designer  for  a  very  much  directed  and 
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purposeful  avenue  to  do  the  redesigns  efficiently  so  that  at  least  one  perhaps 
only  a  few  redesigns  will  finish  the  job. 

The  final  point  that  must  be  noted  is  that  both  types  of  reduction  concepts 
discussed  in  this  section  and  in  Section  4.4,  that  is  the  efficient,  modal  con¬ 
trol  reduction  and  the  efficiency- filtering  can  be  used  simultaneously.  The 
truncation  of  modal  structural  states  based  on  the  criterion  of  this  section  to 
retain  a  2nR-order  structural  dynamics  model  can  be  viewed  as  a  closed-loop  model 
reduction.  On  the  other  hand,  the  truncation  of  efficiency-states  to  obtain  a 
control  input  as  the  feed  back  of  a  reduced  number  of  rip-ef f iciency  states  can 
be  viewed  as  the  closed-loop  controller  reduction.  The  control  law  design  can 
thus,  be  designed  for  a  2nR-th  reduced-order  design  model  in  the  form  of  reduced 
rip- ^-states '  feedback,  thus,  affecting  a  simultaneous  model/controller  reduction 
technique . 

In  the  following  section,  we  discuss  a  number  of  additional  aspects  to 
guide  the  designer  in  using  the  efficiency  approach.  The  process  readily  pro¬ 
duces  more  eificient  designs  as  we  will  present  in  the  next  section. 

4 . 6  Further  Considerations  for  Efficient  Model/Controller  Reductions 

Certain  other  features  of  the  approaches  discussed  in  ths  chapter,  must  be 
noted.  Consider  the  pair  e,  sq  for  the  design  model.  Their  components  ei-cjA® 
and  sq*— cJ(l-Aj)  have  the  property  that  the  indices  i  in  the  set  {e1}  with  the 
highest  rip  components  will  not  necessarily  correspond  to  the  indices  i  in  the  set 
(sq1)  with  the  lowest  np  components.  This  anomaly  depends  on  the  relative  values 
of  the  coefficients  cf  which  are  normalized  square  initial  disturbances  along  the 
principal  controller  directions.  As  this  may  seem  disheartening,  one  cannot, 
however,  change  the  fact  that  efficiency  or  the  inefficiency  of  the  structure 
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control  system  can  bt-  a  strong  function  of  the  initial  disturbance  x0.  Some 
disturbances  are  better  handled  by  a  given  system  design  whereas  others  result 
in  poor  efficiency  performance. 

In  choosing  the  np  efficiency  components  ei(  we  do  not  necessarily  consider 
the  highest  rip  of  them  although  this  seems  like  a  logical  choice.  Further  care 
must  be  taken  to  identify  the  set  {e^  i-1,  .  .  .  .rip) .  A  principal  efficiency  com¬ 
ponent  et  may  be  high  either  because  of  a  high  cf  and  low  A®,  or  low  cf  and  high 
A®  or  both.  Recalling  that  a  spillover  quotient  sqi  is  given  by  (3  ''*5'',  in  the 
first  case,  the  corresponding  spillover  coefficient  sqA  will  also  be  high,  in  the 
second  case,  sqt  will  be  low  and  in  the  third  case  it  will  be  relatively  lower. 
All  depends  on  relative  separations  of  cf  and  Af,  so  retaining  a  high  ei  may  also 
mean  retaining  a  high  spillover  quotient  sqi .  The  word  of  caution  is  that  before 
dominant  et's  are  taken  into  account  corresponding  sq^'s  must  also  be  scru¬ 
tinized.  Similar  caution  must  be  taken  in  case  of  e*  and  sqf. 

Along  the  same  lines,  if  an  initial  control  design  has  a  low  efficiency, 
it  makes  little  sense  to  select  modes  based  on  their  relative  contributions  to 
an  already  low  efficiency.  Instead,  one  should  analyze  their  contributions  to 
a  high  spillover  quotient  and  retain  or  delete  them  based  on  that  account.  Con¬ 
versely,  one  should  bring  in  new  modes  into  the  control  design  model  that  were 
discarded  in  the  original  analysis. 

The  model  reduction  process  that  is  being  discussed  here  seems  to  assume 
that  the  controller  configuration  (number  of  inputs  and  their  locations)  is 
fixed.  In  reality,  however,  model  reduction  (or  mode  selection)  process  and  the 
controller  configuration  are  interdependent.  For  design  purposes,  one  would 
normally  try  to  achieve  a  desired  level  of  efficiency  for  a  reduced  model  order 
of  2nR  and  with  an  upper  bound  on  the  number  of  inputs  m  <  m,^.  Hence,  because 
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the  efficiencies  are  also  a  function  of  controller  configuration,  the  controller 
reduction  process  should  also  utilize  the  controller  configuration  changes.  The 
ultimate  goal  is,  for  the  given  initial  disturbance,  to  identify  the  best  reduced 
order  model  with  the  best  input  configuration  to  go  with  it  to  achieve  a  high 
efficiency  for  the  structure-control  system. 

In  the  above  discussion,  it  is  the  initial  disturbance  represented  in  the 
efficiency  coefficients  cf  responsible  for  not  matching  the  set  (ei)  vs.  (sq^J 
of  (ef)  and  (sq*).  Hence,  for  a  given  control  design,  selection  of  the  best 
principal  controller  directions  will  change  if  the  initial  disturbance  changes. 
An  alternate  approach  to  selecting  the  set  rip  of  principal  controller  directions 
which  avoids  the  initial  disturbance  is  to  consider  the  controller  directions 

associated  with  the  rip  highest  characteristic  efficiencies  {A?  i-1 . rip)  . 

Another  feature  of  this  choice  is  that  the  highest  set  (A®)  will  now  correspond 
to  the  np  lowest  of  characteristic  spillover  values  {Af=2 . np}  . 

Recall  from  Section  3,  that  from  the  properties  of  a  Rayleigh's  quotient 
efficiency  is  bracketed  by  Araax>e>Amin .  It  follows  that  a  control  system  whose 
efficiency  will  be  least  sensitive  to  the  initial  disturbance  x0  must  have  a 
minimal  separation  between  Araax  and  Amin,  that  is  if  the  characteristic  effi¬ 
ciencies  are  dense  in  a  short  interval.  Ideally,  if  Amax=Arain,  the  control  system 
will  be  uniformly  efficient,  better  said  equalized  in  all  principal  directions 
and  will  have  the  same  efficiency  regardless  of  xD.  This  would  be  desirable 
feature.  In  addition,  if  Amin  can  be  made  high  then  the  system  will  also  have 
a  high  efficiency  regardless  of  x0.  The  price  that  must  be  paid  to  make  the 
structure-control  system  efficiency  insensitive  to  x0  is  that  all  principal 
controller  directions  will  now  become  equally  significant,  hence,  truncation  of 
efficiency  states  will  not  be  possible  as  each  £-state  and  its  principal  direc- 
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tion  will  contribute  equally  tc  the  efficiency.  Therefore,  one  will  not  be  able 
to  discriminate  among  the  principal  controller  directions  to  affect  their  trunca¬ 
tion  in  the  £ -space. 

Another  feature  related  to  this  last  observation  is  the  sensitivity  of  the 
stability  of  the  system  to  the  principal  controller  directions.  Since  stability 
of  a  linear  system  does  not  depend  on  the  initial  conditions,  the  stability  of 
a  system  truncated  based  on  the  efficiency  will  depend  not  on  the  magnitudes  of 
c’ne  truncated  efficiency  components  e1(i-np+l,  .  .  .  ,  2n)  but  on  the  modal  properties 
and  tt  associated  with  the  truncated  e^  The  more  significant  the  modal 
properties  of  the  truncated  principal  controller  directions  are,  the  easier  it 
will  be  to  destabilize  the  system  via  efficiency  truncation.  Again,  because  a 
poor  xQ  can  associate  poor  component  efficiencies  et  with  significant  controller 
modes,  efficiency  reduction  based  on  truncation  of  controller  modes  of  poor  et 
will  conflict  with  maintaining  closed-loop  stability.  Stability  robustness, 
then,  to  efficiency-state  truncation  should  be  good  if  the  truncated  controller 
modes  of  poor  et  are  not  significant.  On  the  other  hand,  if  the  insensitivity 
of  the  efficiency  of  the  structure-control  system  to  initial  disturbance  is 
desired,  since  this  will  necessitate  equalization  of  the  characteristic  effi¬ 
ciencies,  each  controller  mode  will  become  equally  important.  Consequently, 
insensitivity  to  xQ  of  the  efficiency  of  the  system,  in  the  presence  of  a  desire 
to  controller  mode  truncation,  will  come  at  the  expense  of  system  with  a  sensi¬ 
tive  stability. 

In  all  of  these  endeavors,  the  efficiency  and  spillover  quotients  and  the 
modal  analyses  serve  as  the  design  directors  qualifying  and  quantifying  the 
merits  of  design  changes  to  achieve  the  best  possible  performance  with  high  effi¬ 
ciency.  In  the  next  section,  a  number  of  illustrations  are  presented  utilizing 
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all  that  is  discussed  in  this  chapter.  Again,  the  particular  stiuctuie  is  the 
ACOSS-4  tetrahedral  structure  used  in  the  previous  chapters. 

4 . 7  Illustrative  Examples  on  ACOSS-4 

The  same  situation  as  in  Section  3.6  is  considered. 

Decomposition  of  Efficiency  Quotients 

For  the  sake  of  brevity,  the  decomposition  of  the  numerator  and  denominator 
control  power  matrices  are  presented  only  for  a  4-th  order  control  design  model. 
The  efficiency  components  et  and  coupled  and  decoupled  component  efficiencies  of 
the  structural  modes  are  also  shown  in  Appendix  C. 

Example  4.1:  This  is  the  same  case  as  in  Examples  3.1.  The  numerator/denomina¬ 
tor  control  power  decomposition  matrices  [N]  and  [ D >  are  given  on 
in  Appendix  C  for  both  global  and  model  efficiencies.  Efficiency 
components  are  eA ,  e*  and  component  efficiencies  ej ,  e*  for  coupled 
('Method  1)  and  decoupled  (Method  2)  cases  are  also  shown. 

Example  4.2:  The  case  of  Example  3.4  is  given  with  the  same  information  listed 
under  Example  4.1. 

Efficiency  Decomposition  Using  nr-principal  Controller  Modes 

Example  4.3:  The  same  system  of  Example  4.2  is  considered.  However,  only  the 
six  (np«6)  highest  characteristic  global  and  relative  model  effi¬ 
ciencies  were  used  in  the  decomposition  of  numerator/denominator 
control  power  matrices  and  the  component  efficiencies  of  the 
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structural  modes  were  based  on  the  6  retained  dominant  controller 


di tec t ions . 


In/Ef f ic iencv  -  f iltering;  Controller  Reduction 

Example  4.4:  Example  3.5  with  n-8  (natural  order)  modes  and  m-4  was  considered. 

The  original  model  had  e*-30.75%,  e-80.65%,  sq*-0.628  and 

sq-19.35%.  Although  the  model  efficiency  e  is  good,  the  global 
efficiency  e*  seems  relatively  poor.  Hence,  an  inefficiency  of 
70%  is  associated  with  the  global  efficiency  whi^h  is  reflected  in 
the  value  of  sq*.  (Note  that  by  definition  sq**l-e*) .  Since 
there  is  more  inefficiency  associated  with  e*  than  efficiency,  the 
modal  analysis  was  done  on  sq*.  The  efficiency- filtering  was 
performed  on  the  modes  of  sq*  by  truncating  the  highest  rip  global 
spillover  efficiency  components  sq*  i-1 , 2 ,  .  .  .  ,rip .  Up  to 
spillover  components  could  be  filtered  without  destabilizing  the 
original  control  design.  The  original  gain  matrix  G  was  simply 
projected  through  G*GT2L2-G(  1  -TjLj^  ,  as  in  Eq.  (4.41)  to  obtain  the 


filtered 

gain 

matrix 

G. 

Note  that  the 

np  principal  global 

spillover 

components 

were 

deleted, 

not 

retained,  so  it 

is 

necessary 

to 

retain 

the 

remaining 

components  represented 

by 

subscript  2.  Conversely,  the  projection  (1-T1L1)  eliminates  the 
highest  contributing  states  represented  by  subscript  1. 
Simulation  results  for  np=4  and  Up-8  are  shown  by  figures  SIM32 
(Figures  4.1,  4.2  and  4.3)  and  SIM34  (Figures  4.4,  4.5  and  4.6), 
respectively.  The  efficiencies  of  both  of  these  sq*- filtered 
cases  were  recomputed.  For  rip-2 ,  the  new  efficiency  quotients 
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were  e*-2f>  /1%  .  e—  +.32%;  for  n^-4  the  new  values  were  e*~'*0 . 9% , 
e«73  18>  Hence,  sq* - f i 1 ter ing  did  improve  the  global  efficiency 
over  the  original  design  by  about  1/3.  Other  pertinent  details 
are  shown  on  the  computer  outputs  listed  in  Appendix  C.  Note  that 
the  simulation  results  show  hardly  any  difference  from  the  origi¬ 
nal  design  SIM10  (Figures  3.10,  3.11,  and  3.12).  The  consequence 
of  this  particular  study  is  that  the  16th  order  design  model  can 
be  controlled  by  the  feedback  of  only  8  efficiency  states,  hence, 
resulting  in  a  50*  reduction  in  the  controller  size. 

Efficient  Model  Selection  (Reduction) 

Example  4.^  The  model  efficiency  of  the  16-th  order  control  design  model  of 
the  first  8  natural  modes  with  four  inputs  (n=*8  natural,  m=4) 
which  was  studied  in  Example  3.5  and  Example  4.4  was  e=80.25%.  In 
an  effort  to  improve  the  relative  model  efficiency  e  further  by 
selecting  a  different  set  of  8  modes  other  than  the  lowest  8 
natural  modes,  the  component  efficiencies  of  Example  3.5  were 
scrutinized.  The  component  efficiencies  for  both  e  and  e*  are 
listed  in  Appendix  C.  Modes  1,  2,  7  and  8  have  the  lowest  coupled 
component  efficiencies  e*  and  ej  according  to  Method  1.  Hence, 
these  modes  were  replaced  with  the  modes  that  were  not  in  the 
original  set.  A  redesign  was  done  for  the  new  16th  order  design 
model  with  the  new  selected  set  n=( Modes:  3-6,  9-12).  The  model 
efficiency  increased  to  e=97.21%,  although  the  global  efficiency 
was  reduced  to  e*-10.27%.  Simulation  results  are  shown  in  Figures 
SIM35  (Figures  4.7,  4.8  and  4.9). 
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Figure  4.5: 


Suboptimal  system  LOS,  n-8  (modes:  1-8),  m-4  (#1-4),  example 
4.4;  controller  reduction  via  efficiency  filtering  for  nP-8 
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Evaluation  model  LOS,  n-8  (modes:  3-6,  9-1?’  *»•- (#1  4), 
example  4.5:  efficient  model  selection  (reduction) 
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5.0  EFFICIENCY  OF  SYSTEMS  WITH  COMPENSATOR  DYNAMICS 


5 . 1  Introduction 

The  developments  presented  in  the  previous  chapters  were  based  on  a  static 
state  feedback  control  law  F-Gx.  In  this  chapter,  we  consider  the  case  whereby 
the  feedback  control  law  is  obtained  via  a  dynamic  compensator  and  extend  the 
efficiency  concepts  studied  in  Chapters  2-4  to  consider  the  effect  of  compensator 
dynamics.  The  extension  of  the  previous  results,  in  this  case,  is  straightfor¬ 
ward  and  the  idea  of  compensator  efficiency  within  the  structure-control  system 
arises  naturally. 

In  Section  5.2,  efficiency  quotients  are  revisited  including  the  compen¬ 
sator  dynamics.  Two  definitions  of  compensator  efficiency  are  given  following 
the  component  efficiency  definitions  introduced  in  Section  4.  A  most  common 
dynamic  compensator  is  the  state-estimator  in  the  control  loop.  In  Section  5.3, 
we,  therefore,  pay  special  attention  to  the  design  of  a  reduced-order  state 
estimator  relevant  to  the  subject  of  this  research.  In  particular,  an  rip-th 
order  efficiency  state  estimator  design  is  given  where  rip  efficiency  states 
are  estimated  to  implement  the  reduced-order  control  law  F^Gfe!  discussed  in  Sec¬ 
tions  4.3  and  4.4.  Here  np  is  the  number  of  dominant  principal  controller 
directions  which  will  have  been  identified  from  an  efficiency  modes  analysis  of 
the  system,  hence,  np<2n.  Finally,  in  Section  5.4  we  consider  formally  the  quan¬ 
titative  and  qualitative  effect  of  control  power  spillover  on  the  response  of  the 
residual  dynamics  not  considered  in  the  2n-th  order  control  design  model. 
Although  this  topic  may  have  been  discussed  in  Chapter  2,  we  postponed  its  pre¬ 
sentation  until  after  the  discussion  of  compensator  dynamics,  so  that  it  will 
seal  appropriately  all  the  developments  presented  heretofore  and  conclude  this 
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investigation  by  underlining  the  critical  nature  of  the  efficiency  of  structure- 
control  systems. 

No  illustrative  examples  are  presented  in  this  ehiptei 


x  -  AcL  x  ,  x  -  [xT  x*]T  (5.5) 

A+BHcC  |  BGC 

Ao,  =  -  (5.6) 

BcHcC+KC  |  AC+BCGC 
»  — 

and  the  control  law 

F  -  Cx  ,  G  =  [HcC  Gc]  (5.7) 


Since  the  definitions  of  the  global  efficiency  and  the  model  efficiency  are 
defined  in  terms  of  only  inputs  F  to  the  dynamic  system  not  on  how  they  are 
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computed  or  derived  their  definitions  Eqs.  (3. 6-3. 8),  are  still  valid,  hence, 
using  the  form  of  F  as  given  by  Eq.  (5.3)  in  Eqs.  (3.6),  (3.7)  we  obtain 

SR  =  J FTRRFdt  -  xJFxq  (5.8) 

SM  -  /FrRgFdt  -  xSPgio  0.9) 

where  P  are  the  (n+nc) -dimensional  power  matrices  obtained  from  the  Lyapunov 
Equations : 

AjPR  +  PrAcL  +  GtRrG  -  0  (5.10) 

AJPM  +  PmAcL  +  GtR£G  -  0  (5.11) 

in  which  Rr  and  R^  are  still  as  given  by  Eqs.  (3.6)  and  (3.7) 
rR  =  DtM~  1D  ,  R^  =  DTEnE^D 

The  model  efficiency  e  and  the  spillover  quotient  sq  are  computed  exactly  as 
described  in  Chapter  3  by  using 

PN  -  PM  ,  PD  -  **  (5.12) 

and  all  that  is  discussed  in  Chapters  3  and  4  will  apply  equally  well  to  control 
with  compensator  dynamics.  The  global  efficiency  e*,  however,  needs  a  little 
careful  consideration.  Recall  that  this  definition  involved  control  of  the  DPS 
with  a  continuously  distributed  input  profile,  where  the  optimal  control  was  the 
independent  feedback  of  the  modal  states  of  the  DPS  structure.  An  important  fea¬ 
ture  of  this  optimal  DPS  control  is  that  it  is  dynamically  similar  to  the 
dynamics  of  the  control  design  model  of  the  structural  modes  in  the  closed-loop 
system.  Hence,  even  in  the  case  of  a  dynamic  compensator,  we  shall  insist  on  the 
similarity  of  the  DPS  distributed  control  design  to  the  structural  modes  alone, 
regardless  of  the  modes  in  the  closed- loop  system  due  to  the  additional  compensa¬ 
tor  dynamics.  Thus,  if  {p}  —  {ps,  pc)  is  the  eigenvalue  spectrum  of  the  closed- 
loop  system,  Eq .  (5.5),  with  the  compensator,  as  a  union  of  the  structural  eigen¬ 
value  spectrum  { ps)  and  the  compensator  eigenvalue  spectrum  { pc ) ,  the  control 
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power  S*  of  the  optimal  distributed  solution  will  be  found  by  using  (ps)  in  the 
equations  given  in  Section  2.4. 

The  global  efficiency  will  then  be 
S* 

e*  _  -  (5.13) 

SR 

Although  at  this  point,  it  is  not  clear  whether  S*  is  greater  or  less  than  5R 
with  a  dynamic  compensator  the  quotient  e*  will  still  preserve  its  function  as 
a  practical  concept  since  it  will  relate  the  power  performance  of  the  actual 
design  represented  via  SR  to  an  ideal  baseline  design  S*  being  a  measuring  stick. 
Specifically,  the  function  of  the  compensator  will  be  highlighted  in  such  a 
setting.  Thus  writing  the  control  power  S* 

2S*  =  xJP*  x0  (5.14) 

the  global  efficiency  e*  will  be  given  by 

XqP*x0  (5.15) 

e* - 

x3PRi0 

where  P*  is  again  computed  as  in  Appendix  B  for  the  spectrum  Ip*)  =  lps). 
Alternately,  augmenting  P*  with  zeros 


e*  can  be  written  as 

XqP*X0 


P*  = 


p*  a. 


Pnrxn  ^n^xn 


e* 


(5.17) 


Next,  partition  the  power  matrices  and  PR  corresponding  to  the  structural 
modal  states  x  and  the  compensator  states  xc 


(5.18) 


where  PR  and  P^  are  the  nxn  upper  left  partitions  of  PR  adn  PM  obtained  from  the 
(nHic) -dimensional  I.vanunov  RniiaMnnc  (5.10  and  5.11).  In  the  case  of  e*  the 
care  that  had  to  be  observed  in  defining  it  by  Eq .  (5.13)  resolves  itself  auto¬ 
matically  regardless  ot  augmentation  of  P*  with  zeros  in  Eq .  (5.16)  and  it  would 
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be  plausible  to  require  independent  feedback  of  structural  modal  states  to  yield 

P*. 

On  the  other  hand,  if  the  compensator  initial  states  are  non-zero,  the 
efficiencies  e  and  e*  will  be  different  than  what  will  be  obtained  by  Eqs . 
(5.20).  Indeed,  a  careful  choice  of  Xqq  can  enhance  the  efficiency  of  the 
control  system  with  the  participation  of  coupling  terms.  The  case  studies  we 
presented  in  earlier  in  Chapter  4  demonstrated  that  coupling  terms  in  the  control 
powers  could  be  beneficial. 

This  brings  us  to  the  point  of  defining  compensator  efficiency  of  the  con¬ 
trol  system  as  a  contribution  to  the  overall  efficiencies  e  and  e*.  This  can  be 
done  easily  considering  the  developments  presented  in  Section  4  on  component 
efficiencies.  The  concepts  presented  therein  can  be  generalized  to  this  case  by 
associating  the  subscript  j  in  the  equations  with  the  structural  modal  states  and 
compensator  states,  hence,  we  consider  j-x,c  and  rewrite  Eq.  (4.1)  as: 

e  -  Eej  -  e8  +  ec  (5.21) 

j 

where  ej-es  and  ej-e,.  follow  from  Eqs.  (4.16)  and  (4.17) 


NDS  +  NSD  NDC  +  NCD 

es  -  -  .  ec -  (5.22) 

2D2  2D2 

Ds  “  Trace [ PRX0 ] s  ,  Dc  -  Trace[PRX0Jc  (5.23) 

Ns  -  Trace  [PMX0]S  ,  Nc  -  Trace  [P'X],.  (5.24) 

“  xoxo  (5.25) 

Similarly  for  e* 
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“  es  +  e 


(5.26) 


e* 


2 

j 


N*D*  +  N'D 


2D2 


N*D„  -f  N*D 


2D2 


(5.2/) 


N*  =  Trace [P*XC]  =  S*  -  N*  ,  N*  -  0 

Ds ,  Dc  and  D  are  as  given  by  Eqs .  (5.23),  (5.24)  and  (5.19).  We  shall  refer  to 
ec  and  e*  as  the  compensator  efficiencies. 

Alternately,  decoupled  compensator  efficiencies  can  be  defined  as  in  Sec¬ 
tion  4.2.  It  is  easy  to  see  that,  in  this  case 


j-s,c  (5.28) 


*oPsxo 

,  nc  -  >9;oPcxco 

XqP§xO 

>  “  xcoPcxco 

x5p*x0 

,  n*  -  0 

(5.29) 

(5.30) 


and  the  expected  efficiencies  are  as  in  Eq .  (4.19) 

ns  +  nc  n* 

E(e)  -  — - -  ,  E(e*)  -  —  (5.31) 

ds  +  dc  ds  +  dc 

which  rely  on  the  expectation  that  the  cross-coupling  terms  between  the  struc¬ 
tural  and  compensator  states  will  average  out  to  zero. 

We  must  emphasize  that  the  developments  resented  here  do  not  depend  on  the 
specific  control  theory  used  cu  design  the  controls.  In  fact,  the  efficiency  of 
the  system  will  serve  to  evaluate  the  control  design  in  a  most  essential  way. 
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For  example,  the  compensator  can  again  be  designed  as  optimal  compensators  in  the 
sense  of  Linear  Quadratic  Gaussian  (LQG)  designs  or  as  fixed  order  optimal 
compensators  via  Minimum  Energy  Optimal  Projection  (MEOP)  technique  (Ref.  14). 
The  efficiency  of  all  such  design  will  still  need  to  be  considered  as  addressed 
in  this  investigation. 

5 . 3  Reduced-Order  Efficiency-State  Estimators 
The  feedback  control  law 
F  =  Gx 

as  given  in  the  previous  chapters  requires  the  knowledge  of  structural  modal 
states  of  the  2n-th  order  control  design  model.  Typically,  a  full  order  state- 
estimator  might  be  used  to  implement  the  feedback  law,  Eq.  (3.10)  by  using  an 
estimate  x  of  the  structural  states.  The  2n-th  order  state-estimator  constitutes 
a  dynamic  compensator  discussed  in  Section  5.2,  Eqs .  (5.2  and  5.3)  and  has  the 
form 


x  -  Ax  +  BF  +  Ky  (5.32) 
F  -  Gx  (5.33) 

where  A  denotes  estimator  matrices  to  be  found  and  G  is  the  mx2n  dimensional  con¬ 
trol  gain  matrix  designed  for  the  system  by  any  suitable  approach.  The  solution 
for  the  estimator  in  terms  of  the  compensator  notation  of  Section  5.2  is  known 
to  be 

A  =  A-KC  -  Ac  ,  B  -  B  -  Bc  ,  G  -  Gc  ,  Hc  -  0  (5.34) 
where  K  is  the  estimator  gain  chosen  to  assign  a  desired  set  of  estimator  (com¬ 
pensator)  eigenvalues  (pc).  In  this  case,  the  separation  principle  holds  and  the 
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4-nth  order  closed-loop  eigenvalues  will  be  { p }  -  (ps  pc)  where  ps  are  the 
eigenvalues  of  the  closed-loop  matrix  Acl-A+BG  for  the  structural  dynamics. 

Hence,  all  of  the  developments  in  Section  5.2  apply  to  this  full -state 
estimator  and  the  impact  of  the  estimator  design  on  the  overall  system  efficiency 
can  be  studied  by  following  the  procedure  given  therein.  In  this  vein,  the 
estimator  efficiency  concept  surfaces  as  the  compensator  efficiency. 

In  practice,  instead  of  full-order  estimator,  a  reduced-order  estimator  of 
order  nc“2n-my  can  be  designed  where  my  is  the  number  of  measurements.  This 
procedure  is  standard  text  book  material  and  does  not  deserve  any  special 
attention  here. 

Another  estimator  which  would  be  especially  relevant  to  the  theme  of  this 
investigation  would  be  an  efficiency- state  estimator,  that  is  an  estimate  t  of 
the  efficiency  state  e  introduced  in  Chapter  3.  Since  the  efficiency  state  e  is 
a  2n- dimensional  vector,  a  full-order  e -state  estimator  will  not  be  of  interest. 
Since  it  will  merely  be  related  to  the  full-order  modal  state  estimator  by  the 
efficiency  modal  matrix 

x  -  Te  (5.35) 

On  the  other  hand,  a  reduced-order  estimator  in  the  e -space  will  be  of  interest 
if  we  recall  the  efficiency-filtering  approach  discussed  in  Section  4.4.  We 
showed  in  Section  4.4  that  a  reduced-order  control  law,  Eq .  (4.35)  F-Gffj  can  be 
used  where  are  the  dominant  efficiency  states  identified  along  the  rip<2n 
principal  controller  directions.  In  Section  4.4,  it  was  also  shown  that  the 
control  law  was  equivalent  to  a  projection  (filtering)  of  the  2n-dimensional 
structural  states,  implying  essentially  a  full-order  state  feedback  of  the  form 
F  -  GTjLjX 
as  given  by  Eq.  (4.38). 
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However,  the  reduced-control  law  can  be  implemented  in  reality  as  the  feed¬ 
back  of  np  retained  efficiency  states  if  an  estimate  of  these  states  are 
available  from  an  np-th  order  efficiency-state  estimator.  It  turns  out  that 
design  of  such  a  reduced  np-th  order  r.-estimator  is  possible.  Recalling  that 
ei  and  x  are  related  via 
£i  -  Lxx 

where  Lx  is  the  npx2n  upper  partition  of  the  T'1  matrix,  we  recognize  that  el  are 
rip  linear  combinations  of  the  structural  modal  states.  An  estimator,  which  can 
estimate  these  specific  linear  combinations  of  the  x-states  will  be  precisely 
what  is  needed  to  implement  the  reduced  order  efficiency  state  feedback  control 
law  (4.35).  The  design  of  estimators  of  linear  combinations  of  states  was 
presented  in  Ref.  (15)  in  a  different  context  and  can  be  adapted  in  our  case  to 
design  an  rip-th  order  ^-estimator . 

Assume  the  -estimator  in  the  form 

«l  -  A  ea  +  BF  +  Ky  (5.36) 

A,  B  and  K  are  npXnp,  npxm  and  npxnty  dimensional  matrices  to  be  found.  Multi¬ 
plying  the  modal  structural  dynamic  equations  on  the  left  by  L: 

Lxx  -  LjAx  +  l^BF  (5.37) 

Subtract  Eq .  (4.37)  from  Eq .  (4.36) 

e  -  Ljx  -  A(e1-L1x)  +  (KC-L1A+AL)x  +  (B-L1)BF  (5.38) 

where  the  term  ALxx  was  added  and  subtracted  on  the  right  hand  side.  If  K  and 
B  are  chosen  so  as  to  satisfy 
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(5.39) 


B  =*=  Lj 

KC  -  I^A  +  AL  -  0  (5. AO) 


and  introducing  the  estimation  error  vector 

£E  =  £i  -  Lxx  (5 .41) 

Eq.  (5.38)  reduces  to 

£E  -  A£e  (5.42) 

which  has  the  solution 

At 

£  e  ”  e  £  eo  .£eo“£e(0)  (5.43) 

Hence,  if  the  estimator  matrix  A  is  chosen,  such  that  it  has  eigenvalues  with 
negative  real  parts,  then  €E(t)=>0  from  which  we  obtain 

t-MO 

lim  -  lim  Lxx 

t-*™  t~>CO 

and  the  will  be  estimated  exactly  as  the  desired  combination  of  states  x  for 
large  t. 

If  we  make  the  simple  choice 


A  -  diagfaq J  aq  <  0  i-l,...,np 


(5.44) 


equation  (5.40)  can  be  solved  for  the  estimator  gain  matrix  K  by  rewriting  it  as 
2n-simultaneous  equations 

CtKt  -  AtL|  -  L|At  (5.45) 
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where  the  2nxmy  CT  matrix  is  assumed  to  have  rank  niy,  ny<2n.  Denoting  the  niy- 
dimensional  row  vectors  of  the  estimator  gain  matrix  K  by  kj  and  the  rip  column 

vectors  of  the  right  hand  side  of  Eq.  (5.45)  by  ai(  i-1,2 . rip,  we  write  Eq. 

(5.45)  in  the  form 

CTki  -  a£  i-1,2, ...,np  (5.46) 

The  solutions  to  Eq. (5.46)  can  be  obtained  by  singular  value  decomposition  of  the 
measurement  matrix  CT 

CT  -  UZV*  (5.47) 

where  U  and  V  are  2nx2n  and  myXitiy  unitary  matrices;  *  denotes  complex  conjugate 
transpose,  and  Z  is  the  2nxmy  matrix  of  singular  values 


Z 


Z0  -  diag  ejj  j-1,2 . niy 


(5.48) 


The  solution  for  the  estimator  gains  k*  can  be  shown  to  be 

kx  -  VZ^USai  i-1 ,  .  .  .  , tip  (5.49) 

where  Uj  is  the  myX2n  upper  partition  of  U*.  It  follows  that  the  estimator  gain 
matrix  K  will  be  given  by 

Kt  -  VZq1  Ui(ATL|-L?AT).  (5.50) 

This  completes  the  design  of  ^-estimator . 


5 . 4  Residual  Dynamic  Response  to  Control  Power  Spillover 

Consider  the  control-design  model  and  the  uncontrolled  system  dynamics 

x  -  Ax  +  BF  F  -  Gx  (5.51) 

Xu  ”  AyXu  +  BuF  (5.52) 

The  response  of  the  closed  loop  control  design  model  is 

x  -  ^cL(t.t0)Xo  ,  *(t)  -  eAcLt  (5.53) 

CL 
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where  <f>CL  is  the  closed-loop  transition  matrix  for  the  control  design  model. 
Denoting  the  control  period  over  which  Xq  is  regulated  by  T=(tf-t0),  the 
spillover  excitation  of  the  uncontrolled  modes  can  be  measured  by  the  total 
square  response  during  T.  The  total  squared  modal  displacement  response  of  the 
uncontrolled  dynamics  will  be  given  by 

Yg  -  S  f[  £2udt  (5.54) 

r=n+l  Lo 

The  uncontrolled  modal  displacements  due  to  control  spillover  excitation  are 
given  by  the  convolution  integrals 

t 

(ru  “  “ribru/  S inwru ( t - r ) Gx ( t ) dr  (5.55) 

to 

Substituting  (5.55)  and  (5.53)  into  (5.54),  we  obtain 

t 

Yu  “  Xo[J  ^cl(t'  .t0)B5Qy(t,r,r')Bu^CL(r,to)drdrdt]xo  (5.56) 

co 

where  B0  is  the  input  influence  matrix  of  the  truncated  dynamics  as  given  in  Eq. 
(2.29).  The  time  dependent  weighting  matrix  Qy  is  given  by 

Qy  -  [w“2]  [ Sz( t , r , t  ' )  ]  (5.57) 

[wuz]  -  diag  [w“z]  r-n+1 . 

[ S2(  t .  r  ,  r '  )  ]  -  diag  [Sinwru(t-r  '  )Sinwru(t-r)  ] 

The  weighting  matrix  Qy  obviously  involves  the  cross  correlations  of  the  impulse 
response  functions  of  the  residual  modes  and  is  inversely  proportional  to  the 
square  of  the  uncontrolled  natural  frequencies.  Equation  (5.56)  will  then, 
noting  that  after  integrations  effectively  an  additional  w'2  will  evolve  so  that 
the  results  will  we  inversely  proportional  to  yield 
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(5.58) 


Y§  -  x£Uy(a>u4)  xo 

where  Uy  must  be  a  positive  definite  resultant  matrix  of  the  weighted  integration 
in  Eq.  (5.56)  since  Yg  is  positive  definite  and  the  argument  u„4  is  shown  to 
highlight  the  nature  of  dependency  of  the  results  on  residual  natural 
frequencies . 

Next,  considering  the  spillover  control  power  S0  defined  by  Eq.  (2.2g)  and 
the  Eqs.  (5.51-5.53) 

Su  *  Xo[/tf/l  ft  ^cl  (T  '  » t0)  Gtb3b0G<^cl  (r,t0)drdr’dt]Xo  (5.59) 

o  o  o  g 

which  must  yield 

Sy  *  xSPtfx0  _  PR  .  pM  (5.60) 

where  the  power  matrices  PR  euJ  ?£  are  available  as  required  for  efficiency 
analysis.  If  the  control  law  is  derived  via  a  dynamic  compensator,  that  is  if 
Eqs.  (5.5)  and  (5.7)  are  used  instead  of  (5.51),  as  described  in  Section  5.2,  it 
is  straight  forward  to  see  that  the  counterparts  of  Eqs.  (5.58)  and  (5.59)  are 
obtained  by  replacing  xc  with  x0,  <f>Ch  with  the  ^CL,  the  closed-loop  transition 
matrix  of  Eq .  (5.50)  and  G  with  G.  The  weighting  matrices  BjjBu  and  BjjQyBu 
will  remain  the  same.  Consequently ,  Pu-PR-Pc  will  take  the  place  of  Pu  in  Eq. 
(5.60). 

Hence,  recognizing  that  the  above  derivation  will  remain  the  same  quantita¬ 
tively  whether  a  compensator  dynamics  is  included  or  not  we  return  to  Eq.  (5.58), 
(5.59)  for  Yy  and  Su-  Comparing  Y0  with  Su  we  note  that  Y0  is  a  weighted  version 
of  the  control  power  wasted  in  the  uncontrolled  modes  where  the  weighting  matrix 
is  Qy  for  Y§  and  unity  for  S0.  Hence,  the  residual  response  and  the  wasted  con- 
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trol  power  are  correlated.  The  significance  of  this  result  is  that  a  measure  of 
residual  dynamic  response  is  available  through  the  wasted  power  Sy,  which  is  in 
turn  available  from  the  efficiency  analysis  directly  without  ever  involving  the 
residual  dynamics  in  the  computations.  It  follows  that  within  the  efficiency 
concepts  presented  in  this  investigation,  evaluation  model  response  simulations 
are  not  necessary  to  infer  degradation  of  the  control  design  model  behavior. 

Since  Qy  is  inversely  proportional  to  the  fourth  power  of  uncontrolled 
natural  frequencies,  it  follows  that  if  natural  frequencies  larger  than  unity  are 
in  the  uncontrolled  dynamic  model  the  net  effect  will  be  that  the  matrix  Uy  in 
Eq .  (5.58)  will  be  considerably  smaller  than  the  control  power  matrix  P{j.  There¬ 
fore,  it  may  take  considerable  control  power  spillover  or  considerable  ineffi¬ 
ciency  in  the  control  design  model  before  any  significant  degradation  of  response 
of  control  design  model  due  residual  excitation  can  be  detected.  This  is  cer¬ 
tainly  not  surprising  as  it  is  well  known  that  high  frequency  modes  are  difficult 
to  excite,  hence,  considerable  control  power  is  bled  by  those  modes  even  for 
minimal  excitation  of  them  if  the  structure  control  design  is  not  efficient.  To 
put  it  more  bluntly,  it  is  more  critical  to  align  the  controller  directions  with 
the  structural  modes  in  the  design  model  to  obtain  high  efficiencies  if  higher 
modes  are  not  in  the  control  design  model. 

The  above  analysis  proves  the  assertion  made  in  Section  2.8  that  perform¬ 
ance  of  the  structural-control  system  must  not  be  judged  on  response  studies 
alone.  Although  no  residual  dynamics  response  can  be  detected  in  the  system  out¬ 
puts,  the  controller  can  still  be  wasting  control  power  operating  inefficiently. 
In  space  missions,  ultimately,  the  power  used  inefficiently  by  the  control  system 
has  to  be  generated  on  board  the  spacecraft  and  therefore,  has  the  potential  to 
impact  the  design  of  power  subsystem  and  other  subsystems  and  even  the  mission 
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significantly.  The  efficiency  results  presented  up  to  this  point,  will  establish 
a  most  important  interdisciplinary  link  of  controller  design  to  other  important 
disciplines  of  space  systems  design.  A  good  structural  control  design  must 
therefore,  strive  for  an  efficient  controller  which  can  be  accomplished  only 
after  an  efficiency  analysis  is  carried  out. 

Furthermore,  if  an  inefficient  controller  is  used  to  control  a  model  with 
high  frequencies  in  the  design  model  which  will  require  considerable  control 
power  and  less  than  unity  low  frequency  inodes  are  in  the  uncontrolled  set,  it 
would  take  only  a  small  amount  of  inefficiency  to  produce  considerable  residual 
mode  excitation  and  degradation  of  the  system  response.  In  such  cases  design  of 
efficient  controllers  will  become  even  a  more  critical  task. 
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6.0  OTHER  APPLICATIONS  OF  CONTROLLER  EFFICIENCY  &  APPLICATION  TO  CSDL-ACOSS-6 


6 . 1  Introduction 

In  this  chapter,  we  present  some  illustrations  of  the  variety  of  structure- 
control  system  designs  that  can  be  addressed  via  the  efficiency  concepts.  Speci¬ 
fically,  in  following  the  subject  of  Section  5.4  on  "Residual  Dynamic  Response 
to  Control  Power  Spillover,"  we  give  a  vivid  illustration  of  the  performance  of 
an  inefficient  controller  in  Section  6.2. 

Another  important  topic  in  the  design  of  structure  -  control  system  is  to 
determine  a  good  input  configuration,  that  is  how  many  actuators  to  use  and  where 
to  put  them,  as  well  as  quantifying  the  effectiveness  of  actuators.  In  Section 
6.3,  we  give  illustrations  of  how  efficiencies  can  be  used  to  select  the  input 
configuration . 

Another  application  of  the  efficiency  concept  can  be  to  choose  the 
weighting  matrices  (or  parameters)  in  the  application  of  the  LQR  theory.  Because 
the  efficiencies  are  defined  regardless  of  the  particular  control  theory  used  in 
the  design  of  the  control  law,  the  opportunity  exists  that  the  LQR  weighting 
matrices  can  be  chosen  to  improve  the  controller  efficiency.  This  is  done  ir. 
Section  6.4. 

As  a  final  analysis  and  design  study  by  using  the  efficiency  technique,  in 
Section  6 . 5  we  present  results  of  efficiency  modes  analysis  of  a  representative 
large  flexible  spacecraft  configuration,  the  Charles  Stark  Draper  laboratory 
AC0SS-6  (Model  2)  structure. 

6 . 2  Inefficient  Control  with  Truncated  Low  Frequency  Modes 

It  was  shown  in  Section  5.4,  that  the  total  square  modal  response  of 
truncated  dynamics  of  a  system  was  a  weighted  version  of  the  control  power  spill- 
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over  and  the  weighting  was  inversely  proportional  to  the  fourth  power  of  the 
natural  frequencies  of  the  truncated  modes.  Hence,  the  efficiency  or  ineffi¬ 
ciency  of  a  controller  becomes  critical  with  low  frequency  truncated  modes,  also 
from  a  response  point  of  view. 

As  an  illustration  of  this,  we  refer  the  reader  to  Examples  3.1  and  3.3 
with  respective  simulations  given  in  figures  SIM39  (Figures  3.1,  3.2  and  3.3)  and 
SIM41  (Figures  3.4,  3.5  and  3.6),  in  Chapter  3.  Both  of  these  controllers  had 
low  efficiencies,  e*=0.173%,  e=4.31%  and  e*=6.2%,  e=28.59%.  The  control  design 
model  was  4-th  order  and  consisted  of  high  frequency  modes  (#11  and  12),  the 
lowest  frequencies  were  left  in  the  truncated  dynamics.  Note  that  the  degrada¬ 
tion  in  the  response  will  be  as  serious  as  the  magnitude  of  the  lowest  frequency 
in  the  truncated  dynamics.  For  the  ACOSS-4  structure  the  lowest  natural  fre¬ 
quency  is  wj-1.34  rad/sec.,  so  it  too  has  an  attenuating  effect  on  the  control 
power  spillover  (proportional  to  1/wJ)  ,  nevertheless,  not  as  strong  as  the  other 
truncated  modes.  Certainly  if  wx<l ,  the  control  power  spillover  effect  would  be 
magnified  in  the  residual  response.  However,  Wj-1 . 34  was  still  low  enough  to 
accentuate  the  residual  response  degradation  due  to  inefficient  controller. 
In  both  SIM39  (Figures  3. 1-3. 3)  and  SIM41  (Figures  3. 4-3. 6)  evaluation  model 
responses  are  significantly  degraded. 

Two  other  examples  of  evaluation  model  response  degradation  due  to  ineffi¬ 
cient  controllers  are  presented  below: 

Example  6.1:  n-4  {Modes:  9-12),  m-2(#l,#2).  The  system  has  e*-8.39%,  e-45.67%. 

SIM28  shows  the  responses  (Figures  6.1  through  6.3). 
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Example  6.2:  n=4  {Modes:  9-12),  m-1  (#5).  The  system  has  e*=0.202%,  e-21.33%. 

SIM30  shows  the  responses  (Figures  6.4  through  6.6).  Clearly  the 
more  inefficient  system  of  Example  6.2  has  serious  degradation  of 
response 

6 . 3  Actuator  Selection  Via  Efficiencies 

The  full-order  12  mode  ACOSS-4  evaluation  model  was  considered  in  conjunc¬ 
tion  with  a  full  set  of  12  actuators,  n=12,  m=*=12 .  Because  there  is  no  model 
truncation,  the  relative  model  efficiency  of  the  designs  are  e-100% .  However, 
e*  is  not  necessarily  100%  since  the  12  inputs  do  not  represent  a  continuously 
distributed  input  field  although  they  may  come  close  to  representing  it.  The 
number  of  inputs  were  decreased  from  m-12  to  m=l ,  in  each  case  e*  was  computed, 
the  changes  in  e*  with  the  deletion  of  each  input  being  a  reflection  of  the 
effectiveness  of  that  input.  The  following  results  were  obtained: 


ACTUATOR  EFFECTIVENESS  ON  ACOSS-4  EVALUATION  MODEL 


m-=Actuator  Number 

%e*Global  efficiencies 

Ae* 

Actuator  Rank 

12 

70.88 

-0.59 

12 

11 

71.47 

0.28 

11 

10 

71.19 

7.13 

3 

9 

64.06 

6.74 

4 

8 

57.32 

13.85 

2 

7 

43.47 

31.73 

1 

6 

11.74 

0.7 

10 

5 

11.04 

0.78 

8 

4 

10.26 

0.72 

9 

3 

9.54 

1.16 

7 

2 

8.38 

4.10 

6 

1 

4.28 

4.28 

5 

First,  note  that  e*  is  valid  indicator  of  the  quality  of  the  input  configuration 
and  the  effectiveness  of  an  actuator  can  be  measured  by  its  contribution  to  e* . 
Hence,  each  actuator  was  ranked  according  to  its  effectiveness.  The  most  effec- 
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tive  5  actuators  were  (#7,8,10,9,1).  For  this  model,  the  above  ranking  of  the 
actuators  is  in  agreement  with  the  actuator  ranking  by  computing  the  work  done 
bv  each  actuator  during  the  control  period  (Ref.  16).  A  second  study  was  done 
by  using  a  reduced-order  model  n=l-8.  In  this  case,  since  e  is  less  than  100%, 
it  too  was  considered  in  judging  the  actuator  effectiveness  along  with  the  global 
efficiency* . 


ACTUATOR  EFFECTIVENESS  ON  ACOSS-4  WITH  A  REDUCED-ORDER  MODEL 


m 

% 

e* 

% 

e 

Rank* 

Ae*% 

Ae% 

Rank 

i 

5  . 

69 

56 

.  92 

8 

5. 

.69 

56. 

.92 

1 

2 

8  . 

41 

55 

.26 

12 

2. 

.72 

-1. 

.66 

10 

3 

25  . 

79 

75 

.67 

1 

17. 

.38 

20. 

41 

2 

4 

30. 

76 

80 

.65 

9 

4. 

97 

4. 

98 

5 

5 

45. 

30 

70 

.62 

2 

14. 

.  54 

-10. 

03 

12 

6 

53. 

09 

60 

.84 

5 

7. 

.79 

-9. 

78 

11 

7 

60. 

48 

76 

.66 

6 

7. 

.39 

15. 

82 

3 

8 

67. 

i_n 

00 

84 

.21 

7 

7. 

.37 

7. 

,55 

4 

9 

72. 

.40 

83 

.47 

10 

4 

.55 

-0. 

.74 

9 

10 

81 

.19 

87 

.29 

3 

8. 

.79 

3, 

.82 

6 

11 

89 

.  36 

87 

.58 

4 

8 

.17 

0. 

.29 

7 

12 

93 

.08 

87 

.28 

11 

3. 

.72 

-0. 

.30 

8 

Again,  note  that  e*  is  a  consistent  indicator  of  actuator  effectiveness  whereas 
e  is  anomalous,  both  rankings  are  different  and  in  this  case  because  the  control 
design  model  changed,  ranking  of  actuators  according  to  e*  has  also  changed  from 
the  previous  example  set  when  the  evaluation  model  was  considered.  Simulation 
results  of  n=8 ,  m=l ;  and  n=8 ,  m=12  are  shown  in  figures  labeled  SIM8  and  SIM18 
(Figures  6.7  through  6.12). 
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6 . 4  Effect  of  LQR  Weighting  Matrices  on  Efficiencies 

The  global  and  model  efficiencies  were  computed  for  different  state  and 
control  weighting  parameters  q  and  r  as  discussed  in  Section  2.8,  for  the 
weighting  matrices  Wx=block- diag  [ q ) ,  Wc=r[l]  for  a  given  control  design  model 
and  input  configuration,  n=*8  (natural)  and  m=l .  The  results  are  shown  below: 

EFFECT  OF  LQR  WEIGHTING  PARAMETERS 


r 

q 

e*% 

e% 

1.0 

4.1 

5.69 

56.93 

0.1 

4,1 

12.12 

56.93 

°.05 

4.i 

16.89 

56.93 

0.01 

4,i 

27.29 

56.93 

1 

i 

4.99 

56.93 

1 

0.5 

4.99 

56.93 

1 

0.1 

4.99 

56.93 

1 

10 

5.26 

56.93 

1 

100 

7.81 

56.93 

1 

1000 

31.86 

56.93 

SIM20  and  SIM26  (Figures  6.13  through  6.18)  show  the  responses  for  r=0.05  and 
q-100,  respectively.  The  above  example  shows  that  the  model  efficiency  e  is 
insensitive  to  control  and  state  ghting,  which  might  be  expected  since  e  is 
the  ratio  of  two  control  powers  SM,  SR  for  the  same  control  design,  hence, 
changing  q  and  r  changes  both  by  the  same  proportion.  On  the  other  hand,  e*  is 
based  on  two  dynamically  similar  control  designs  but  with  different  input 
profiles  hence,  e*  should  be  sensitive  to  changes  in  q  and  r. 

These  examples  are  based  on  a  rather  simple  variation  of  the  LQR  weighting 
matrices.  The  effect  of  changes  in  the  full  weighting  matrices,  for  example,  in 
the  case  of  different  frequency  shaped  weighting  matrices,  may  produce  not  so 
.simple  trends.  In  the  above  examples,  the  insensitivity  of  the  model  efficiency 
coefficient  to  weighting  matrices  can  be  used  to  advantage  in  the  sense  that  it 
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will  be  possible  to  adjust  the  response  time  of  the  system  by  changing  q  and  r 
without  changing  (a  good)  model  efficiency. 

6 • 5  Application  to  a  Large  Space  Structure:  CSDL-ACOSS- 6  (Model  2) 


Next,  a  realistic  large  space  structure  was  considered.  The  structure 
shown  in  Fig.  (6.19)  is  the  Charles  Stark  Draper  Laboratory  model  known  as  ACOSS- 
6  (Model  2)  (ref.  17).  The  finite  Element  evaluation  model  that  was  considered 
was  produced  by  NASTRAN  at  the  Flight  Dynamics  Laboratory  at  WRDC/WPAFB.  The 
model  that  was  studied  for  efficiency  had  294  degrees  of  freedom.  Efficiency  of 
two  different  control  design  models  employing  different  sets  of  inputs  were 
considered.  There  were  available  21  actuators  located  at  positions  shown  in 
reference  18.  Efficiency  of  two  designs  were  analyzed. 

Example  6.3:  In  this  case,  the  control  design  model  had  8  modes,  n=(12,  13,  17, 
21,  22,  24,  28,  30)  with  4  inputs  m~(18,  19,  20,  21).  Actuator 
numbers  are  in  consistence  with  the  designation  listed  in  Ref. 
(18).  The  modes  chosen  were  designated  significant  tor  LOS 
studies  (ref.  19).  The  simulation  results  of  square  x-y  plane 
deflection  at  node  37  are  shown  in  figures  labeled  SIMV1  (Figures 
6.20  and  6.21).  The  global  efficiency  was  e*=0.176%  and  the 
relative  model  efficiency  was  e-0.657%. 


Example  6.4:  The  control  design  model  had  10  modes  n={7-16)  and  10  inputs  m-{l- 
10).  The  same  deflection  as  in  example  6.3  is  shown  in  figures 
labeled  SIMV2  (Figures  6.22  and  6.23).  The  efficiencies  of  this 
system  were  e*=0.236%  and  e=0.401%. 
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Because  Che  first  six  modes  obtained  from  the  structural  eigenvalue  problem  were 
not  reliable  none  of  them  were  taken  in  the  control  design  models  in  both  of  the 
examples.  Such  poorly  known  modes  have  no  consequences  for  the  reliability  of 
efficiency  analysis  as  long  as  they  are  in  the  truncated  dynamics  of  the  system 
as  discussed  in  Section  2.6.  The  computer  results  of  efficiency  modes  analysis 
are  given  in  Appendix  C. 

In  examples  6.3  and  6. A  for  the  ACOSS-6  structure  the  responses  of  the 
evaluation  model  were  not  shown  since  this  would  require  simulation  of  a  588  th 
order  state  model,  clearly  beyond  any  reasonable  expectation.  However,  in  light 
of  the  illustrations  and  theoretical  developments  presented  thus  far,  the  con¬ 
trols  designed  for  the  ACOSS-6  structure  are  extremely  inefficient.  Almost  99.5% 
of  control  power  will  excite  the  truncated  dynamics.  To  make  this  ineffieicney 
even  more  critical,  there  are  some  modes  with  natural  frequencies  less  than  unity 
left  in  the  truncated  dynamics.  Hence,  as  demonstrated  in  Example  6.2  for  the 
ACOSS-4  structure,  in  this  case  for  ACOSS-6  too,  serious  degradation  of  the 
responses  will  occur  in  the  actual  (evaluation  model)  structure  instead  of  the 
control  design  model  responses  shown  in  Figures  (6.21)  and  6.23).  However,  no 
simulation  of  the  evaluation  model  is  necessary  based  on  the  qualitative  con¬ 
clusion  that  the  efficiency  approach  affords.  (It  is  understood,  however,  that 
a  simulatable  model  larger  than  the  control  design  model  and  smaller  than  the 
588-th  order  evaluation  model  can  be  used  to  underline  these  statements.  The 
challenge  here  is  to  resist  this  temptation). 

The  efficiency  obtained  above  for  a  realistic  large  space  structure  are 
extremely  low  indeed  showing  the  need  and  demonstrating  the  merits  of  a  full- 
scale  assessment  of  the  need  for  an  efficiency  analysis.  As  shown  in  Appendix 
A,  the  efficiencies  defined  in  this  research  are  literally  the  power  efficiencies 
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of  the  control  system.  A  0.6%  model  efficiency  simply  means  6  energy  units/sec. 
out  of  1000  were  consumed  by  the  controller,  the  rest  simply  being  wasted.  The 
wasted  power  will  have  been  derived  on  board  the  spacecraft  to  serve  no  useful 
control  purpose.  Such  inefficiencies  in  the  system  cannot  be  tolerated.  It 
should  be  recalled  that  although  the  control  design  model  is  a  much  smaller 
model,  the  information  that  is  obtained  from  them  in  the  form  of  power  efficiency 
pertains  to  the  behavior  of  the  full -order  evaluation  model. 
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Figure  6.2: 


Suboptimal  system  LOS ,  (modes:  9-12)  m—2  (#1,2),  exampl 

6.1;  response  of  an  inefficent  system  e*-8.39%,  e-45.67% 
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Figure  6.3: 


Evaluation  for  Model  LOS,  n-4  {modes:  9-12),  m-2  (#1-2),  example 
o.I;  response  of  an  inefficient  system  e*~8.39%,  e-45.67% 
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Figure  6.6: 


Evaluation  model  LOS,  n-4  (modes:  9-12)  ra-1  (#5),  example  6.2; 
response  of  an  Inefficient  system,  e*-0.2%,  e-21.33% 


00 

co 

CM 

O 

o 

o 

o 

O 

o 

o 

o 

jojjg  J°  0un 


Figure  6.7: 


Globally  optimal  system  LOS,  n-8  (modes:  1-8)  m-1  (#1)  for 
actuator  selection  via  efficiencies 


Time  (sec) 


jcujg  J°  9un 


Figure  6.11: 


Suboptlmal  system  LOS,  n-8  (modes:  1-8),  m-12;  for  actuator 
selection  via  efficiencies 
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Figure  6.12: 


Evaluation  model  LOS,  n-8  (modes:  1-8),  m-12;  for  actuator 
selection  via  efficiencies 
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Figure  6.16:  Globally  optimal  system  LOS,  n-8  (modes:  1-8);  effect  of  LQR 

weighting  on  efficiency,  q-100,  e*-7.81%,  e-56.93% 
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Evaluation  model  LOS,  n-8  (mode.r:  1-b);  effect  of  LQR  weighting 
on  efficiency,  q-100,  e*-7.8U,  e-56.9.1% 


Figure  6.19:  CSDL  -  ACOSS  -  6  (model  2)  structure 
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Figure  6.20:  Globally  optimal  system,  square  x-y  plane  deflection  at  node  37, 

n-8 ,  (modes:  12 , 13 , 17 , 21 , 22 , 24 , 28 , 30} , example  6 . 3 
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Figure  6.22:  Globally  optimal  system,  square  x-y  plane  deflection  at  node  37, 

n-10,  (modes:  7-16),  example  6.4 
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Figure  6.23:  Suboptimal  system,  square  x-y  plane  deflection  at  node  37, 

{modes:  7-16),  m-10  (#1-10),  example  6.4 
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7.0  CONCLUSIONS 


The  concept  of  efficiency  of  an  engineering  system  is  a  time  proven  design 
and  analysis  tool  in  many  disciplines.  S true ture - Control  Systems  should  not  pose 
an  exception  to  this  practice.  This  research  undertook  the  task  of  defining  an 
efficiency  concept  for  structure-control  systems,  and  established  physical 
grounds  for  the  efficiencies  which  would  have  the  potential  to  address  the  issues 
in  the  control  of  large  structural  dynamic  systems.  We  have  shown  that  the 
global  and  relative  model  efficiencies  proposed  and  developed  in  this  research 
are  the  power  efficiencies  of  the  control  system  and  therefore,  they  must  be 
considered  as  fundamental  figures  of  merit  of  a  control  system  design  along  with 
other  more  familiar  performance  measures.  No  less  important  a  feature  of  the 
efficiencies  is  that  the  performance  of  the  evaluation  model  can  be  judged  on  the 
basis  of  che  control  design  model  alone  without  explicitly  involving  the  full- 
order  evaluation  model,  a  feat  made  possible  only  by  a  unique  judicious  choice 
of  the  weighting  matrices  involved  in  the  efficiency  definitions. 

The  investigation  developed  the  new  concept  of  controller  modes  signifying 
characteristic  avenues  of  efficient  control  power  utilization.  These  efficiency 
modes  are  complementary  to  the  familiar  concept  of  structural  modes  and  together 
they  determine  the  effectiveness  of  any  structural  control  design.  Ultimately 
what  counts  is  the  degree  of  matching  of  these  two  sets  of  modes.  The  efficiency 
modes  analysis  leading  to  the  recognition  of  principal  controller  directions  led 
to  a  procedure  to  give  structure  to  the  internal  workings  of  the  structural 
control  system  through  the  control  power  matrices.  We  have  attempted  to  exploit 
this  internal  characterization  of  the  system  by  defining  component  efficiencies 
and  efficiency  components.  Furthermore,  we  proposed  and  developed  model/con¬ 
troller  reduction  approaches  based  on  that  structural  information  by  correlating 
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the  structural  m.d<  a  anu  the  controller  (efficiency)  modes.  We  have  illustrated 
the  significance  of  the  efficiency  concepts  by  a  host  of  examples  on  realistic 
structures.  The  end  point  of  this  investigation  offers  many  other  avenues  of 
research.  To  mention  only  two:  optimization  of  structural  control  system  for 
efficiency  subjected  to  other  design  constraints,  design  of  control  gains  by  a 
direct  approach  utilizing  the  system  efficiency  as  the  target  design  criterion. 
We  have  demonstrated  that  the  Linear  Quadratic  Regulator  Theory  can  produce 
highly  inefficient  control  designs  although  they  are  regarded  as  optimal. 
However,  the  optimality  of  an  LQR  design  is  subjective  and  theory  merely  serves 
as  a  means  of  obtaining  control  gains  in  some  sense.  On  the  other  hand,  the 
definitions  of  power  efficiencies  of  the  control  system  constitute  absolute  non- 
dimemsional  unique  performance  measures.  No  subjectivity  is  involved  in  their 
definitions.  The  efficiency  of  all  control  systems  regardless  of  the  methodology 
by  which  they  are  designed  can  be  analyzed  by  the  approaches  developed  in  this 
research  and  their  merits  can  be  established  on  a  common  ground. 
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APPENDIX  A:  PROOF  OF  FUNCTIONAL  SR  AS  CONTROL  POWER 

The  control  evaluation  functionals  SR,  SM  and  S*  which  are  uniquely  defined 
for  the  efficiency  analysis  of  structure-control  systems  represent  power  quan¬ 
tities,  true  rate  of  change  of  energy,  a  fundamental  concept  in  the  design  of 
engineering  systems.  To  see  this,  consider  the  starting  point  of  this  investiga¬ 
tion  in  Chapter  2  which  is  the  control  functional  SR  defined  over  the  DPS: 

SR  -  J  fT(p,t)m_1(p)f(p,t)dD(p)dt  (A— 1) 

corresponding  to  the  partial  differential  equations  of  motion. 

m(p)u(p,t)  +  L(up , t)  -  f(p,t)  (A-2) 

First,  perform  a  dimensional  analysis.  In  Eq.  (A-2)  m(p)  has  the  units  of 
mass/dimD  where  dimD  is  the  dimension  of  Domain  of  the  Structure,  mass  is  the 
mass  or  inertia  distribution  over  D.  For  one-dimensional  beam  elements 
dimD=length,  for  plate  elements  dimD  is  an  area  etc...  Consistently,  the  input 
function  f(p,t)  represents  a  force  density  over  the  domain  D,  it  has  the  dimen¬ 
sion  Load/dimD  where  load  can  be  a  force  or  a  moment.  Substituting  the  dimen¬ 
sional  equivalents  of  the  terms  appearing  in  the  definition  (A-l) 

Load  mass 

f(p.t)  -  - —  ,  m(p)  =  - 

dimD  dimD 

Load  dimD  Load 

SR  =  J  - —  — —  -  dimD  *  d(time) 

dimD  mass  dimD 

Load2  *  Time 

SR  - - - -  (A— 3) 

mass 

and  noting 
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displacement 

Load  -  mass  *  acceleration  -  mass  *  - 

time2 


Expression  (A-3)  yields  after  simplification 
Load  *  displacement  Energy 

SR  - -  -  -  Power  (A— 4) 

Time  Time 

Since  SM,  S*  are  also  derived  from  SR  by  an  invariant  frame  transformation,  SM 
and  S*  are  also  power  quantities. 

Next,  we  consider  the  mechanical  aspects  of  this  power.  To  make  this 
easier  to  see,  consider  the  one— dimensional  mechanical  analogy  of  the  DPS 
equation  of  motion  (A-2) ,  the  familiar  spring  mass  system  where  the  spatial 
dimension  has  been  (integrated)  out. 

k 

mu  +  ky  -  f(t)  ,  -  -  w2  -  rad/sec2  (A-5) 

m 

for  which  SR  is 

SR  -  Jra-1f2(t)dt  (A-6) 

Multiply  Eq.  (A-5)  by  m_1f(t)  and  integrate  over  time,  and  recognize  that  the 
right  hand  side  becomes  SR: 

Jfudt  +  Jw2fudt  -  SR  (A-7) 

substituting  u  =  du/dt 
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Jfdu  +  J"cj2fudt  —  SR 


(A-8) 


Again  the  left  hand  side  has  power  units  and  specifically  the  first  integral  on 
the  left  is  the  area  under  the  load  versus  velocity  diagram.  Alternately,  write 
the  equation  of  motion  (A-5) 

m(u+cj2u)  —  f(t)  (A— 9) 

and  define  the  effective  rigid  body  acceleration  that  the  forcing  function  f(t) 
imparts  to  the  mass  m  as  uR  =  u  +  w2u  where  uR  is  the  effective  rigid  body  dis¬ 
placement  due  to  f(t).  Hence,  (A— 5)  is  equivalent  to 
muR  =  f(t) 

Multiply  this  last  expression  by  uR 

muRuR  =  f(t)uR  (A-10) 

d 

H - mu2  =  fuR  (A— 11) 

dt 

from  which  we  recognize  the  effective  rigid  body  kinetic  energy 
TR  =  hmujj 

and  Eq .  (A-10)  is  the  more  familiar  power  expression 

d 

- Tr  -  fuR  -  P 

dt 

The  incremental  power  change  associate  with  Eq .  (A— 12)  is 

fdup  =  dP 

The  effective  rigid  body  incremental  velocity  duR  can  be  written  as 
dup  =  (u  +  o;2u)dt 


(A— 12) 


(A-13) 
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Multiplying  on  the  left  by  f  and  recognizing  (A-12))  we  obtain 
Jf(u+a>2u)dt  -  JfduR  -  JdP  -  P 

where  P  is  the  total  mechanical  power.  The  integral  on  the  left  is  SR  by  virtue 
of  Eq .  (A-8) ,  hence,  we  get 

SR  =  Jm_1f2dt  =  JfduR  <=  P 

We  have  thus  also  shown  mechanically  that  SR  functional  is  a  total  power  expres¬ 
sion  associated  with  the  control  input  f  in  trying  to  move  the  structure  with  an 
effective  rigid  body  velocity  UR.  It  must  be  recognized  that  uR  is  not  a  fic¬ 
titious  mathematical  expression.  Both  uR  and  uR  as  derivatives  are  well  defined 
quantities  at  all  times  in  terms  of  u  as  given  by  Eq.  (A-13) .  The  concepts  of 
effective  rigid  body  velocity/acceleration  are  also  plausible  because  this  is 
what  the  controller  attempts  to  impart  to  the  structural  mass  via  the  thrust  of 
control  input  f.  The  controller  actuators  at  any  instant  will  not  recognize  that 
the  mass  has  flexibility.  The  effect  of  flexibility  is  inherent  in  the  structure 
defined  by  the  domain  of  m,  and  it  manifests  itself  to  result  in  an  actual  accel¬ 
eration  u  by  extracting  w2u  from  uR 

uR  -  oj2u  “  u  (A-14) 

which  provides  the  mechanism  to  bleed  off  control  power  into  the  flexible 
dynamics  to  store  it  as  elastic  energy  in  the  structure. 

Next,  we  consider  some  specific  control  hardware. 

Thruster  Jets 

For  a  thruster  using  an  exhaust  jet,  the  thrust  force  f(t)  is  the  control 
input  and  given  by 

f  -  Ap(t)  Vex  (A— 15) 
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where  Vex  is  the  exhaust  jet  velocity  and  m(t)  is  rate  of  flow  of  propellant  mass 
nip.  Form  the  control  functional  SR  for  such  a  thruster 

SR  =  J'm“1f2(t)dt  (A-161 

and  substituting  (A-15)  into  (A-16)  by  assuming  constant  exhaust  velocity 

V2 

ex 

SR  =  f  m|(t)dt  (A— 17) 

m 

in  which  the  propellant  flow  rate  is  the  control  variable.  For  simplicity,  the 
change  in  the  Lotal  mass  in  due  to  propellant  flow  was  neglected.  The  kinetic 
energy  rate  of  flow  carried  away  by  the  exhaust  jet  is  on  the  other  hand 

d 

£p  -  H - KV2J  =  HV2xi?,p(t)  =  P(t)  (A-18) 

dt 

where  P  is  the  power  of  the  exhaust  jet  at  time  t.  Substituting  from  (A-18)  into 
(A-17)  for  iiip,  we  get 

4 

SR  -  -  JP2(t)dt  (A-19) 

mV«x 

Denote  the  root  mean  square  power  in  the  exhaust  beam  by  Prms  and  write 
PLeT  -  f P2(t)dt 

where  T  is  the  control  (regulation)  time.  Next,  introduce  an  average  character¬ 
istic  power  Pc  for  the  combination  of  the  structure  to  be  controlled  and  the 
nature  of  the  thruster  propellant 

m&^sp 

Pc  *  hmV2x/T  =  - - -  (A-20) 

2T 
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where  Isp  is  the  specific  impulse  of  the  thruster  propellant  and  gc  is  the  gravi¬ 
tational  acceleration  on  earth.  Now  the  SR  can  be  rewritten  as 


:r 


sK  -  v 


>i\. 


41  k\ 


SK  -  t 


)P, 


Wogolap 


VA-  2i) 


in  which  WD  is  the  on  earth  weight  of  the  structure. 

Equations  (A-21)  clearly  shows  that  SR  is  proportional  to  the  root  mean 
square  power  of  the  exhaust  beam  ivwer  with  a  nondimens ional  proportionality 
constant  which  provides  an  interesting  combination  of  control,  structure  and 
thruster  propellant  parameters. 


Proof-Mass  Actuators 

As  another  controller  example,  we  consider  a  proof-mass  actuator  of  the 
type  discussed  in  Ref.  (20).  The  actuator  dynamics  given  in  Ref.  (20)  accounts 
for  coupling  with  the  structural  dynamics  and  typically  represents  a  compensator 
dynamics,  a  subject  which  was  discussed  in  Section  5.2.  For  simplicity  of  this 
discussion,  we  shall  disregard  the  interaction  of  the  actuator  with  the  structure 
and  represent  the  force  f(t)  it  exerts  on  the  structure  by  the  equation 

f(t)  -  KampKfV(t)  (A-22) 

where  V(t)  is  the  control  voltage,  K£  is  a  force  coefficient  (force/amper) .  The 
Ramp  is  given  in  terms  of  the  armature  current  and  voltage: 

1 

Wa)  -  Iain<8)/V(s)  =  -  (A-23) 

Ls  +  R 

s  is  the  LaPlace  variable,  L  and  R  are  the  inductance  and  resistance  of  the 
actuator.  If  the  time  constant  L/R  is  small,  then  Kamp-  R"1.  In  this  case  the 
actuator  force  output  is  simply  given  by 
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(A-24) 


f(t)  =  KfIarm(t'  =  KffTMt) 

Forming  SR 

S8-  Jm  1 12(.  t  ) ctt  -  Jm  lKjR  2V2dt  -  Jin  >K2l2dt 

Multiplying  and  dividing  by  R  and  recognizing  that  12R  is  the  electrical  power 

P  consumed  by  the  actuator  to  control  the  structure  we  write, 

K2  K2T 

S8  -  - jRI2dt  =  - Pm 

mR  mR 

for  which  we  defined  the  mean  power  Pm=  T_1jRI2dt  representing  the  average 
electrical  power  of  the  proof-mass  actuator  over  the  control  period.  The 
coefficient  of  Pn  can  also  be  shown  to  be  a  nondimens ional  quantity  relating  the 
control,  structure  and  actuator  parameters  as  in  the  case  of  thrust  jets 
discussed  earlier. 

The  quantities  SM  and  S*  by  definition  have  the  same  units  as  SR  and 
therefore,  the  efficiency  quotients  defined  e,  e*  in  this  report  are  physically 
the  power  efficiencies  of  the  controller. 
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APPENDIX  B:  CLOSED  FORM  SOLUTION  FOR  S* 


The  S*  is  the  control  power  for  the  continuously  distributed  input  field 
over  DPS  where  the  modal  inputs  are  based  in  the  independent  modal— space  control 
approach  (Regs.  1,  4).  The  modal  state-space  dynamics  for  the  2n-th  order  design 
model  is: 


x 


0 

1 


fr*(t) 


r=l , 2 , . . . ,n 


( B— 1 ) 


with  the  modal  states  xr=[£l£r]T,  and  f*(t)  is  the  modal  input  given  by 


"  glut  +  Sr  2^  r 


Sri  ^r  Pz\Pt2  >  Sr2  **  ( P  tlPzl^ 


prl  and  pc2  are  the  r-th  closed-loop  modal  eigenvalues  corresponding  to  the  open- 
loop  structural  eigenvalues  ±jwr. 

The  control  power  S*  as  defined  in  Chapter  2  can  now  be  written  as 

s*  -  Jf;2(t)dt  =  s  xjoPX 

r  *1 


P'  is  the  2x2  optimal  control  power  obtained  as  the  solution  of  n  2x2  independent 
Lyapunov  equations 


A?p;  +  p;ac  +  g;tc;  =  0  r-l . n  (B-2) 

0  l] 

Ar  -  !  •  g;  =  [gr\  gr*2] 

H  °J 

The  solution  of  these  modal  Lyapunov  equations  is  easily  obtained  as 
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Pi  po 


P5  P’z 


CP— 3) 


wsig| 


2(gi-w2) 


g2(gi“w2) 


“  glg2 


(B-4) 


(B-5) 


2(gi-w2) 


2g2(gi“w2) 


(B-6) 


in  which  g^g^  ,  g-^z'tZ  >  r-l,2,...,n. 
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APPENDIX  C:  COMPUTER  RESULTS  OF  EXAMPLES  ILLUSTRATED  IN  SECTIONS  3,  U  &  6 


Nomenclature  for  the  computer  outputs 

The  following  consists  of  description  oi  terminology  encountered  in  the 
computer  program  that,  may  not  be  self  explanatory.  Some  output  statements  will 
not  be  of  any  interest  to  the  reader  as  they  are  only  for  program  control  and 
debugging  purposes.  Hence,  no  explanation  is  given  for  them. 

*  Reduced  principal  efficiency  components  NP  is  the  number  of  effi¬ 
ciency  states  considered,  rip-1 , 2  ,  .  .  , NP . 

*  Structural  natural  frequencies  for  the  evaluation  model  are  ordered 
in  increasing  order.  The  first  n  are  in  the  control  design  model  in 
increasing  order,  the  rest  are  in  the  truncated  model  in  increasing 
order.  NE  is  the  evaluation  model  order.  N  is  the  controlled  struc¬ 
tural  modes. 

*  Init.  Gen.  Displacement  and  Init.  Gen.  Velocity  are  q(t0)  and  <5(t0) 
for  the  evaluation  model. 

*  S*  is  S*  ,  SM  is  S”  and  SR  is  SR. 

*  Principal  Components  Analysis  of  _  _  efficiency  is  _  efficiency 

modes  analysis. 
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PO  is  the  dimensional  efficiency  states  initial  conditions  e  .  CE 
is  squared  efficiency  coefficients  cf.  i-=l ,  2  ,  .  .  ,  NP  ,  Eq  .  (3.36). 

Characteristic  efficiencies  are  Af  in  random  order.  Efficiency  vec¬ 
tors  are  the  efficiency  components  ei ,  total  efficiency  is  the  effi¬ 
ciency  (not  in  percent). 

Eff .  .  .  Nuin.  .  .  Control  Cost  is  XqPnx0 
Eff . . .  Dcnora. . .  Control  Cost  is  XqPdx0 

Numerator/Denominator  Component  cost  is  (XqP^)^  j-l,2,..,n  for  each 
structural  mode. 

Unordered ...  Comp  Eff  (Num.  comp.  Costs/Denom)  is  (X0PN)jj/' XoPDx0 . 

Unordered  Comp ,. Eff iciencies  is  Equation  (4.16). 

Method  1  is  the  "coupled  efficiencies"  cf.  Equation  (4.16)  of 
structural  degrees  of  freedom. 

Method  2  is  the  "decoupled  efficiencies"  of  Equation  (4.17)  of 
structural  degrees  of  freedom. 

Mode  (ir.  modsel)  gives  in  decreasing  order  coupled  efficiencies  e^  for 
structural  modes.  In  modsel  implies  that  the  modes  are  desigr^ted  as 
labeled  in  selected  modes  for  control  design.  For  example  if  14th 
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and  31st  modes  are  selected  for  control  in  modsel  they  are  designated 
as  the  first  and  2nd  selected  modes,  respectively. 

*  NPth  order  efficiency  (ENP)  indicates  that  efficiencies  are  computed 
via  Eq .  (4.20)  -  (4.22)  by  considering  NP  efficiency  components. 

*  Block-diag  numerator  and  denominator  Comp,  costs  are  side  by  side  n^ 
and  dj  in  Eq .  (4.17).  Each  row  corresponds  to  a  j-th  structural  mode. 

*  Block-diag.  eff .  .  .Num/Denom.  .  .Control  cost  is  the  numerator/denomina¬ 
tor  of  Eq .  (4.19).  Unordered  comp,  efficiencies  is  Eq .  (4.17)  for 

each  j  . 

*  Block-diag .. Comp .. Eff . .  is  e^  in  Eq .  (4.17)  in  decreasing  order  for 
the  modes  in  selected  modes. 

*  NPth  order  block-diag.  efficiency  ENP  is  Eq .  (4.19). 
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